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SECTION  I 
INTRODUCTION 

Accelerations  of  less  than  one  second  duration  occur 
in  a variety  of  aerospace  environments  including  escape  system 
and  aircraft  and  surface  vehicle  crash  environments.  The  6570th 
Aerospace  Medical  Research  Laboratory  (AMRL) , Bionics  and 
Biodynamics  Division,  Impact  Branch  (BBI)  has  been  engaged  in  the 
conduct  of  experimental  tests  to  assess  the  effectiveness  of  a 
variety  of  restraint  systems  and  protective  devices  by  evaluating 
biomechanical  responses  of  manikins,  primates,  and  humans  to 
exposure  to  laboratory  simulations  of  these  environments.  The 
biomechanical  responses  to  such  exposure  were  evaluated  from 
strategically  located  and  oriented  monitoring  devices,  such  as 
accelerometers  and  load  cells  whose  outputs  were  recorded  on 
magnetic  tape,  and  high-speed  motion  picture  cameras  mounted 
off board  and  onboard  the  test  vehicles. 

The  reduction  of  electronically  recorded  data  was 
accomplished  in  routine  fashion  using  accepted  methods  to  derive 
the  Gadd  Severity  Index  (SI)  and/or  Head  Injury  Criterion  (HIC) 
from  triaxial  acceleration  time-histories.  While  these  are 
recognized  indicators  of  the  probability  of  injury  resulting 
from  exposure  to  an  impact  environment,  they  do  not  describe 
the  motion  of  body  segments  other  than  at  the  point,  or  points 
monitored.  To  adequately  describe  the  responsive  motion  of 
the  many  body  segments,  quantitative  analysis  of  time  displace- 
ment data  of  the  anthropometric  points  defining  the  various  body 
segments  was  required.  The  photometric  system  designed  and 
implemented  to  satisfy  this  requirement  was  an  integrated 
data  acquisition,  data  reduction  system. 

Data  were  recorded  on  at  least  three  16mm  cameras 
mounted  offboard,  or  in  a combination  of  offboard  and  onboard, 
depending  on  the  environment  being  simulated.  If  vehicle 
accelerations  were  not  to  exceed  a peak  acceleration  of  20g,  two 
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cameras  were  mounted  onboard  the  vehicle  with  two  others  mounted 
offboard.  If  peak  vehicle  accelerations  in  excess  of  20g  were 
anticipated,  all  cameras  were  mounted  offboard. 

When  the  photo  instrumentation  plan  specified  that  all 
cameras  be  offboard,  one  or  two  primary  data  cameras  (depending 
on  the  vehicle  displacement  during  the  event) , were  mounted 
normal  to  the  plane  of  principal  predicted  motion.  Two  other 
cameras  were  mounted  such  that  their  optical  axes  were  mutually 
perpendicular  to  the  optical  axis  of  the  primary  data  camera. 

Wnen  the  photo  instrumentation  plan  required  the  onboard/ 
offboard  combination,  the  two  onboard  cameras  were  located  and 
oriented  with  respect  to  a vehicle  coordinate  system.  The 
locations  and  orientations  of  these  cameras  varied  from  one  test 
program  to  another  depending  on  the  size  of  the  test  subjects, 
type  of  predicted  motion,  etc.,  with  principal  consideration 
being  given  to  the  fact  that  all  tracked  points  had  to  be  observed 
by  both  cameras  throughout  the  entire  event.  Of  the  offboard 
cameras  one  was  located  with  its  optical  axis  normal  to  the  plane 
of  principal  motion  and  the  other  was  located  such  that  its  optical 
axis  was  aligned  parallel,  or  as  nearly  parallel  as  practical, 
to  the  axis  of  vehicle  motion. 

Prior  to  each  test,  fiducials,  one-inch  diameter  circles 
with  alternating  black  and  yellow  quadrants,  were  applied  in 
accordance  with  the  recommended  practice,  SAE  J138  (see  Figures 
B1  and  B3) . After  the  subject  was  adjusted  to  the  initial 
position,  the  locations  of  the  applied  fiducials  were  measured 
and  recorded.  Reference  fiducials,  installed  on  the  test  fixture 
and  surveyed  at  initiation  of  a test  series,  were  checked  prior 
to  each  test  to  determine  if  any  required  replacement. 

The  data  reduction  phase  of  the  system  involved  digitizing 
the  film  plane  coordinates  of  the  reference  fiducials  and 
tracking  fiducials,  time  synchronizing  of  film  frames, 
calculation  of  conversion  factors,  and  electronic  data 
processing. 
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The  projected  film  frame  coordinates  were  digitized, 
frame  by  frame,  using  a Producers'  Service  Corporation  model  PVR 
film  analyzer  interfaced  to  a type  35  teletype  terminal  (TTY) . 

The  TTY  generated  a listing  and  punched  paper  tape  of  formatted 
digital  values. 

The  digitized  data  were  transmitted  from  the  TTY  tape 
reader  to  the  computer  via  voice  quality  telephone  line.  The 
data  file  was  edited  to  correct  format  or  characters  as  required 
then  saved  on  disk  file  and  copied  to  punched  cards  which  served 
as  the  permanent  file. 

The  system  contained  three  processing  programs.  HIFPD 
was  used  to  process  planar  motion  of  points  as  observed  by 
a camera  mounted  normal  to  the  plane  of  motion.  This  program 
was  initially  written  to  yield  time  histories  of  coordinate  points 
relative  to  the  test  vehicle  and  of  the  test  vehicle  relative  to 
the  range.  It  was  shortly  thereafter  modified  to  derive  relative 
linear  and  angular  velocities  and  accelerations  of  the  points 
and  to  prepare  plots  of  these  data. 

Program  SLED  was  designed  to  take  data  digitized  from 
two  camera  views,  synchronize  the  data  timewise,  and  solve  for 
the  most  likely  points  of  intercept  of  up  to  four  pairs  of  rays 
in  a three-dimensional  rectilinear  coordinate  system.  The 
successful  operation  of  this  program  required  precise  determination 
of  the  coordinates  of  the  focal  point  of  each  camera  and  the 
azimuth,  elevation,  and  roll  angles  of  each  of  the  cameras. 

Physical  measurements  of  these  parameters  proved  inadequate,  thus 
the  third  program  in  the  system,  POOCH,  was  developed. 

Program  POOCH  is  used  to  derive  the  location  and 
orientation  of  a camera  and  a factor  including  the  focal  length 
of  the  camera  lens  and  the  magnification  of  the  projector.  It 
optimizes  these  parameters  based  upon  the  surveyed  coordinates 
of  up  to  20  reference  points. 
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The  cameras  were  operated  at  a nomir.al  film  speed  of 
500  frames  per  second,  with  exposure  time  of  0.4  msec  to  0.8  msec 
depending  upon  predicted  maximum  velocities  of  tracked  points 
and  illumination  intensity.  Film  timing  and  camera  synchronization 
were  provided  for  by  the  use  of  pulsed  LED's,  driven  by  a common 
pulse  generator  at  a rate  of  100  pulses  per  second,  recorded  on 
one  side  of  the  film  just  beyond  the  edge  of  the  image  frame. 

Due  to  the  location  of  the  LED  in  the  film  path,  any  given  pulse 
was  recorded  11  to  13  frames  behind  the  coincident  image  frame, 
depending  on  the  size  of  the  film  slack  loop  between  the  supply 
reel  and  the  film  gate.  To  minimize  the  error  due  to  this 
variance,  a synchronizing  pulse,  recorded  on  the  electronic  data 
recording  system,  was  used  to  fire  a strobe  light,  the  flash  of 
which  was  observed  by  all  cameras. 


* SECTION  2 

ANALYSIS  OF  DATA  FROM  ONE  CAMERA 


Horizontal  Impact  Facility  Photometric  Data  Analysis 
Program  (HIFPD)  is  a digital  computer  program  developed  to  analyze 
the  Hyge  Impact  Facility  Photometric  data  for  Impact  Branch  of  the 
Biodynamics  and  Bionics  Division  of  the  6570th  AMRL.  The  program 
was  compiled  and  executed  on  the  CDC  computers  at  Wright-Patterson 
Air  Force  Base.  The  standard  CALCOMP  plot  package  is  used  to 
plot  data  and  thus  must  be  attached  to  load  and  execute  the 
program. 

This  program  inputs  the  code  sheet  data  and  program  con- 
trol parameters  described  in  Appendix  A and  a maximum  of  300 
(MAXN)  frames  of  X,  Z position  data  for  the  range,  sled,  hip,  knee, 
shoulder,  elbow,  head  point  1 and  head  point  2 for  I TYPE  = 0 
or  range,  sled,  head  point  1 and  head  point  2 for  I TYPE  = 1.  The 
data  card  format  are  also  described  in  Appendix  A. 

The  program  computes  the  following  four  types  of  data 
as  requested  by  the  program  control  parameters : 

(a)  The  input  data  versus  frame  number  and  the  frame  to 
frame  differences  are  printed  in  counts.  The  range  difference  is 
subtracted  from  the  frame  to  frame  differences  for  each  of  the 
seven  parameters.  The  only  value  of  this  difference  data  would 
be  to  spot  errors  in  the  data. 

(b)  The  displacements  (X  and  Z)  of  the  hip,  knee,  shoulder, 
^lbow,  head  point  1 and  head  point  2 relative  to  the  sled  are 
computed,  and  a moving  eleven  point  (NP  = 11)  quadratic  least 
square  fit  is  used  to  smooth  the  data.  These  data  are  also  plotted, 
if  requested  on  the  test  setup  card. 

(c)  The  angles  in  radians  between  the  shoulder  and  hip 
and  between  the  head  point  1 and  head  point  2 are  computed  using 
the  above  smoothed  data.  The  angular  velocity  is  computed  in 
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radians  per  second  using  a moving  11  point  quadratic  fit  of  the 
angle  versus  time  data  (computes  derivative  of  least  squares 
equation) . The  angular  acceleration  is  computed  using  a moving 
eleven  point  quadratic  fit  of  the  velocity  versus  time  data. 

These  data  are  also  plotted  as  requested  on  the  test  setup  card. 

(d)  The  linear  velocity  and  acceleration  data  for  any 
combination  of  the  eight  variables  are  computed  as  requested  on 
the  test  setup  card.  For  example,  the  linear  velocity  and 
acceleration  of  the  head  pt  1 relative  to  the  range,  sled  relative 

to  the  range  or  the  head  pt  1 relative  to  the  sled  can  all  be 

computed.  Note  that  range  relative  to  some  other  parameter  cannot 
be  computed.  To  compute  these  linear  velocity  and  acceleration 
data,  the  X and  Z displacements  are  computed  for  the  variable  of 
interest  relative  to  the  reference  variable.  A moving  eleven  point 
(NP  = 11)  quadratic  least  square  smoothing  function  is  applied 
to  both  the  X and  Z time  histories.  The  resultant  displacement 
in  feet  is  computed  from  this  smoothed  data.  A moving  eleven 

point  quadratic  fit  is  applied  to  this  resultant  data  to  obtain 

the  velocity  in  feet  per  second  and  a second  eleven  point  quad- 
ratic fit  is  applied  to  the  velocity  data  to  obtain  the  accelera- 
tion data  in  feet  per  second  squared  and  in  G’s.  These  data  are 
printed  and  plotted  as  requested  on  the  test  setup  card. 

The  three  external  files  used  by  this  program  are  the 
input  file  (unit  5)  used  to  read  all  code  sheet  and  data  cards, 
the  output  file  (unit  6)  used  to  print  all  output,  and  TAPE7 
(unit  7)  used  to  generate  the  plotter  tape.  A magnetic  tape 
must  be  requested  with  TAPE7  as  the  local  file  name. 

The  following  sections  of  this  report  present  a general 
description  of  the  main  proqram  and  all  subroutines  except  the 
CALCOMP  plot  routines.  Flow  charts  are  also  included  for  each 
routine.  Appendix  C contains  a complete  listing  of  the  program 
source  deck  and  Appendix  D contains  a sample  run  complete  with 
all  input  and  output  data  (including  CALCOMP  plots) . 
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2.1 


PROGRAM  HIFPD 


£-  This  main  routine  controls  all  input,  output,  and 

I computations  requested  by  the  test  setup  card  parameters.  All 

| subroutines  required  to  smooth  the  data,  compute  derivatives, 

f and  plot  results  are  called  by  this  routine.  All  program 

j diagnostics  resulting  from  errors  in  setup  or  data  card  formats 

are  printed  by  this  routine. 

Method 

The  program  reads  the  code  sheet  control  cards  described 
in  Appendix  A and  initializes  the  program  print  and  plot  control 
parameters.  The  program  reads  the  card  code,  frame  number,  and 
X and  Z axis  data  for  four  ( I TYPE  =1)  or  eight  (ITYPE  = 0) 
variables  (index  J)  for  each  frame  (index  I)  in  the  test.  Frame 
number  and  card  code  are  checked  for  input  errors;  errors  in 
input  cause  diagnostics  to  be  printed  and  the  test  to  be 
terminated.  If  more  than  MAXN  frames  are  read,  diagnostics  are 
printed  and  all  frames  beyond  MAXN  are  omitted  from  the  analysis. 
The  T ( I ) time  data  are  computed  from  the  frame  number  as  follows; 

T ( I ) = IFR ( I ) /DT 

where  IFR (I)  is  the  frame  number  and  DT  is  the  number  of  frames 
per  second.  If  setup  card  parameter  IRX  is  greater  than  zero, 
the  sign  of  all  X axis  data  are  changed.  After  all  data  are 
read,  a summary  page  is  printed  listing  all  types  of  analyses 
to  be  computed,  printed,  and  plotted  for  this  test. 

When  program  control  parameter  IPR  > 0,  all  raw  input  X 
and  Z axis  data  are  printed  in  counts.  The  frame  to  frame 
difference  data  are  computed  and  printed  for  all  J variables 
from  frames  1 equal  2 to  N as  follows: 

XD ( 1 ) = X(I,1)  - X(I-1,  1) 

XD ( J)  = X(I,J)  - X(I-1,J)  - XD ( 1 ) . 

12 
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XD ( 1 ) is  the  range  difference  from  the  i n frame  and  XD(J)  is 
the  variable  minus  range  difference  for  the  Jth  variable  and  the 
Ith  frame.  The  above  are  also  computed  and  printed  for  the  Z 
axis  data. 

All  X and  Z axis  data  are  adiusted  for  shifts  in  the 
range  reference  reading  and  then  converted  from  counts  to  feet: 


HI  = X(I,1)  - X (1 , 1 ) 

H2  = Z(I,1)  - Z(l,l) 

X(I,J)  = (X (I , J)  - HI)  * CAL(J) 
Z(I,J)  = (Z(I,J)  - H2)  * CAL(J) 


where  CAL(J)  is  the  calibration  factor  for  the  Jtn  variable 
(J  = 2 to  8)  . 

When  program  control  parameters  IPC  < 2 or  IPA  < 2 , 

X and  Z axis  motion  relative  to  the  sled  are  computed  for  variables 
3 to  8 ( or  7 and  8 for  I TYPE  = 1) : 

XD ( I ) = X(I,J)  - X (I , 2) 

ZD ( I ) = Z(I,J)  - Z (1 ,2)  . 

Subroutine  SM  is  called  to  compute  a moving  eleven  point  (NP  = 11) 
quadratic  least  square  fit  to  smooth  the  X and  Z axis  data. 

The  smoothed  data  are  stored  in  arrays  XX(I,JJ)  and  ZZ(IfJJ)  where 
JJ  = J-2.  As  a result  of  the  eleven  point  smoothing,  five  frames 
are  lost  at  the  beginning  and  end  of  the  test  data;  this  is  true 
each  time  the  data  are  smoothed  by  subroutine  SM  or  derivaties 
are  computed  by  subroutine  DERIVl . If  parameter  IPC  < 2,  these 
smoothed  data  relative  to  the  sled  are  printed;  if  IPC  < 1, 
subroutine  CPLT  is  called  to  generate  a CALCOMP  plot  of  X versus 
Z for  all  variables  (J  = 3 to  8) . 
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The  angle  between  the  shoulder  and  the  hip  is  computed 
for  each  frame  using  the  above  smooth  data  when  program  control 
parameter  IPA  < 2.  The  angle  in  radians  is  computed  as 
follows : 

HI  = ZZ  (1 , 3)  - ZZ(I,1) 

H2  = XX (1,3)  - XX (1,1) 

XD (I ) = arctan  (H1/H2) 

where  index  3 is  shoulder  data  and  index  1 is  hip  data  in  the 
XX  and  ZZ  arrays.  Angles  XD(I)  are  adjusted  by  factors  of  2n 
to  make  them  continuous.  Subroutine  DERIV1  is  called  to 
compute  the  angular  velocity  in  radians  per  second  from  a 
moving  eleven  point  (NP  = 11)  quadratic  fit  of  the  XD(I)  data 
and  angular  acceleration  in  radians  per  second  squared  from  an 
eleven  point  quadratic  fit  of  the  velocity  data.  The  angular 
data  are  printed  and,  for  IPA  = 0,  subroutine  CPLT  is  called  to 
generate  CALCOMP  plots  of  the  angular  velocity  and  acceleration 
versus  time  (IP  = 2) . All  above  angular  data  are  computed  in 
a similar  manner  for  head  point  1 minus  head  point  2 data 
(indices  5 and  6 in  arrays  XX  and  ZZ) . 

Parameter  M contains  the  number  of  sets  of  linear 
velocity  and  acceleration  data  to  be  computed  for  one  variable 
(array  ID)  relative  to  another  (array  IR) . For  example,  if 
ID ( 1 ) = 3,  and  IR(1)  =2,  then  for  set  M = 1 the  hip  motion 
relative  to  the  sled  is  computed  for  all  available  frames. 

If  M < 0 and  IPL  < 2,  all  data  for  variables  J = 2 to 
8 are  adjusted  by  subtracting  the  initial  value  as  follows: 

X(I,J)  = X(I,J)  - X(1,J) 

Z ( I , J)  = Z(I,J)  - Z(1,J) 

where  all  X and  Z data  have  previously  been  converted  from 
counts  to  feet.  For  each  of  the  M sets  the  following  are 
computed : 
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JD  = 

ID  (K) 

JR  = 

IR(K) 

DI(I) 

= X (I , JD) 

- X (I , JR) 

DC(I) 

= Z ( I , JD) 

- Z ( I , JR) 

where  K is  the  set  index  (K  = 1 to  M) . When  JR  is  1 above, 

X ( I , JR)  and  Z(I,JR)  data  are  deleted,  because  motion  relative 
to  the  range  has  already  been  removed  from  the  data  (see  calibra- 
tion equations) . Subroutine  SM  is  called  to  apply  the  eleven 
point  quadratic  smoothing  function  to  the  DI(I)  and  DC (I) 
displacement  data;  the  -smoothed  data  are  stored  in  arrays  XD(I) 
and  ZD  (I).  The  resultant  displacement  RES (I)  is  computed 
as  follows: 

RES  (I)  = ^D(I) 2 = ZDd)^. 

The  maximum  XD(I),  ZD (I)  and  resultant  and  the  minimum  XD(I)  and 
ZD (I)  are  also  determined  from  the  above  data.  Subroutine  DERIV1 
is  called  to  compute  the  velocity,  VEL(I) , from  the  resultant 
data  and  the  acceleration,  ACC(I),  from  the  resultant  data  using 
eleven  point  quadratic  smoothing  functions  in  each  case.  The 
acceleration  data  are  converted  from  feet  per  second  squared  to 
G’s: 


ACCG(I)  = ACC (I) /32 . 2 . 

All  the  above  displacement,  velocity  and  acceleration  data  are 
printed  versus  time.  Also,  if  IPL  = 0,  subroutine  CPLT  is  called 
(IP  - 3)  to  plot  the  velocity  in  feet  per  second  and  the 
acceleration  in  G's  versus  time  in  seconds. 

After  all  tests  are  analyzed,  the  CALCOMP  plot  subroutine 
PLOTE  is  called  to  terminate  the  CALCOMP  plot  routines. 
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A flow  chart  for  this  routine  is  shown  in  Figure  1 and  a 
sample  test  complete  with  all  input  and  output  is  listed  in 
Appendix  D. 


Error  Diagnostics: 
Subroutines  Required: 


COMMON  Required ; 


Program  Length; 
Labeled  Common  Length; 
Blank  Common  Length: 


CM  Required  to  Load 
Program: 


YES 


CPLT , DERIV1 , SM,  PLOTS, and  PLOTE 
(PLOTS  and  PLOTE  are  part  of  the 
CALCOMP  plot  package . ) 


The  COMMON  variables  actually  used 
to  transfer  data  are  described  in 
subroutine  CPLT.  Many  of  the  blank 
COMMON  variables  in  this  routine 
are  in  COMMON  only  to  reduce 
required  compilation  storage . 


5342 


8 


23, 


8 

12616 


8 


Approximately  50000 
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2.1.1  Subroutine  CPLT  (T,Y,  Z,  IP) 

Subroutine  CPLT  generates  a CALCOMP  plot  of:  (a)  X versus 
Z displacement  with  respect  to  the  sled  for  all  available  para- 
meters, (b)  time  versus  angular  velocity  and  acceleration,  or 
(c)  time  versus  linear  velocity  and  acceleration  depending  on 
the  value  of  parameter  IP.  The  standard  CALCOMP  plot  package 
is  required  to  load  and  execute  the  program. 

Method 

For  parameter  IP  = 1,  CPLT  generates  one  composite  plot 
of  X versus  Z displacement  for  variable  motion  with  respect  to 
the  sled.  All  X and  Z scaling  are  set  by  CPLT  independent  of  the 
data  and  are  defined  as  follows. 


HIFPD  Flow  Chart 


(a)  The  minimum  X value:  XMIN  = -1.4  - 2.2  * IRX. 

(b)  The  minimum  Z value:  MIN  = 0.0. 

(C)  The  X and  Z increments  per  inch:  DX  = DZ  = 0.4. 

(d)  The  X and  Z axis  lengths:  SX  = SY  = 10. 

The  variables  are  defined  and  thus  plotted  (all  on  one  graph) 
from  indices  II  to  12  when  flag  ICAL(J)  > 0 for  J = 3 to  8. 

Data  to  be  plotted  are  stored  in  blank  common  arrays  XX  and  ZZ. 
All  data  are  checked  to  be  sure  they  are  within  the  plot  scale 
values;  if  not,  they  are  set  equal  to  the  minimum  or  maximum 
value.  Subroutines  LINE  and  SYMBOL  are  called  to  plot  the  data 
and  print  the  legend  on  the  graphs. 

For  parameter  IP  = 2,  CPLT  generates  one  plot  of  time 
(T  or  X array)  in  seconds  versus  angular  velocity  (Y)  in  radians 
per  second  and  acceleration  (Z)  in  radians  per  second  squared. 
The  time  scaling  is  determined  by  the  program  as  follows: 

(a)  the  minimum  time  value  (XMIN)  is  determined  from  the 
initial  time  value,  X(l) , adjusted  to  the  nearest  0.01  less  than 
X(l)  ; 

(b)  the  time  increment  per  inch,  DX  = 0.02; 

(c)  The  time  axis  length  (SX)  is  determined  from  DX  and 
the  total  range  X(N)  - XMIN 

SX  = FLOAT ( IFIX ( (X(N)  - XMIN)/DX)  + 1). 

The  angular  velocity  and  acceleration  minimum  and  increment  per 
inch  scaling  are  set  up  by  calling  subroutine  SCALE  which  checks 
the  data  and  sets  values  accordingly.  The  velocity  scale  is 
printed  on  the  left  side  of  the  graph  and  the  acceleration  scale 
on  the  right  side.  Subroutines  LINE  and  SYMBOL  are  called  to 
plot  the  data  and  print  the  legend  on  the  graph. 

For  parameter  IP  = 3,  CPLT  generates  one  plot  of  time 
(T  or  X array)  in  seconds  versus  linear  velocity  (Y)  in  feet  per 
second  and  acceleration  (Z)  in  G's.  The  time  scaling  is 
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computed  as  per  IP  = 2 above.  The  velocity  and  acceleration  are 
plotted  using  the  same  ordinate  scale.  The  ordinate  length  SX  is 
always  10  inches.  The  minimum  value,  YMIN , and  the  increment  per 
inch,  DY,  are  determined  from  the  data.  DY  will  always  be  10, 
20,  or  30  depending  on  the  total  range  required.  If  the  total 
range  is  greater  than  300,  some  of  the  data  will  be  lost.  All 
data  are  checked  to  be  sure  they  are  within  the  scale  values;  if 
not,  they  are  set  equal  to  the  maximum  or  minimum  value.  Sub- 
routine LINE  is  called  to  plot  the  data  and  subroutines  SYMBOL 
and  NUMBER  are  called  to  print  the  legend  on  the  graph. 

A flow  chart  for  this  routine  is  shown  in  Figure  2 and 
sample  graphs  are  shown  in  Appendix  D. 


Error  Diagnostics:  NONE 

Subroutines  Required;  AXIS,  LINE,  NUMBER,  PLOT,  SCALE, 

and  SYMBOL  (all  are  in  the  standard 
CALCOMP  plot  package) . 


Argument  List:  T - array  of  time  data  points  in 

velocity  and  acceleration  plots. 

Y - array  of  velocity  data  points. 

Z - array  of  acceleration  data  points. 

IP  -flag  which  determines  the  type 
of  graph 

1 - composite  plot  of  variable 

displacement  with  respect 
to  the  sled 

2 - angular  velocity  and 

acceleration 

3 - linear  velocity  and 

acceleration 

Blank  and  Labeled 

COMMON  Variables;  JD  - integer  used  to  print  parameter 

name  in  legend 

JR  - integer  used  to  print  parameter 
name  in  legend 

N - number  of  data  points  used  in 

velocity  and  acceleration  plots 

NP  - number  of  points  used  in  least 
square  fit 
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Setup  absicca 
(time)  parameters 
for  velocity  and 
acceleration  plots; 
call  AXIS  to  draw 
absicca. 


Is 

IP  = 2?' 


Yes 


Determine  maximum 
and  minimum  ordi- 
nate scaling  for 
•’  ear  velocity 
acceleration 
j;  call  AXIS 
| t ' draw  prdinate . 


Check  for  velocity 
and  acceleration 
data  outside  plot 
range. 


Call  LINE  to  plot 
velocity  and 
acceleration  versus 
time  data. 


Call  SYMBOL  and 
NUMBER  to  print 
ID  information. 


Setup  plot  of  X 
versus  Z for  para- 
meter motion  with 
respect  to  the 
sled. 


Setup  ordinate  an 
absicca  scaling  and 
call  AXIS  to  draw 
both  axis. 


Call  PLOTE 
to  shift  plot 
origin  for 
next  plot. 


RETURN 


Figure  2.  CPLT  Flow  Chart. 
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11  - first  point  used  in  composite 

plot 

12  - last  point  used  in  composite 

plot 

XX  - array  of  X axis  displacement  data 

ZZ  - array  of  Z axis  displacement  data 

ICAL  - flag  array  which  identifies 
defined  data 

ICAL ( J)  = 0 - Jth  variable 
undefined 

ICAL ( J)  = 1 - Jth  variable 
is  defined 

HEADL  - array  containing  variable 
names  used  in  legend 

TEST  - test  identification  used 
in  legend 

IRX  - flag  used  to  setup  composite 
plot  X axis  scale 

Subroutine  Length:  1363g 

Labeled  Common  Length:  23g 

Blank  Common  Length;  3456g 

2 . 1 . 2 Subroutine  XM(X,  Y,  YC,  N,  NP) 

Subroutine  SM  is  a smoothing  routine  which  computes 
a quadratic  least  square  fit  of  NP  dependent  variable  data  points 
(Y)  to  compute  each  smoothed  data  point  (YC) . Since  NP  data  points 
are  used  to  compute  each  smoothed  point,  M data  points  are  lost 
at  the  beginning  and  end  of  array  YC,  where 

M = (NP  - l)/2. 


Method 

The  first  (MM)  and  last  (NN)  array  indices  for  which 
YC ( I ) are  computed  are  determined  as  follows: 

MM  = M + 1 
NN  = N - M 
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where  M is  defined  above  and  N is  the  number  of  original  displace- 
ment points  in  array  Y.  Subroutine  QLSQ  is  called  to  compute 
the  , C2 , and  coefficients  for  each  of  the  I smoothed  points 
which  are  then  computed  as  follows: 

YC ( I ) = CL  * X(I)2  + C2*X(I)  + C3. 

A flow  chart  for  this  routine  is  shown  in  Figure  3. 

Error  Diagnostics : NONE 

Subroutines  Required:  QLSQ 

Argument  List:  X = array  of  independent  variable 

Y = array  of  dependent  variable 

YC  = array  of  smoothed  dependent 
variable  data 

N = number  of  original  displacement 
versus  time  data  points 

NP  = number  of  points  used  to  com- 
pute each  smoothed  data  point 

Subroutine  Length:  113g 


2.1.3  Subroutine  DERIV1  (X,  YP,  N,  NP,  ID) 


Subroutine  DERIV1  computes  the  derivative  (YP)  of  the 
dependent  variable  Y.  A quadratic  least  square  fit  of  NP  points 
is  used  to  compute  each  derivative  point;  thus  K points  are  lost 
at  the  beginning  and  end  of  array  UP: 

Where  K = M + M * ID, 

M = (NP  - l)/2. 


ID  = 1 for  first  derivative,  and 
ID  = 2 for  second  derivative. 


Note  that  for  ID  = 1,  array  Y contains  displacement  data  which 
have  already  been  smoothed  using  a quadratic  least  square  fit 
over  NP  points;  thus,  M points  have  already  been  lost  from  the 
original  displacement  data.  For  ID  = 2,  array  Y contains  first 
derivative  (velocity)  data  which  starts  at  array  location 
Y (2*M  + 1) . 
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Method 


The  first  (MM)  and  last  (NN)  array  indices  for  which  YP(I) 
are  computed  are  determined  as  follows : 

MM  = K + 1 
NN  = N - K 

where  K and  M are  defined  above  and  N is  the  number  of  original 
displacement  data  points.  Subroutine  QLSQ  is  called  to  compute  the 
C^,  , and  coefficients  for  each  of  the  I derivative  points. 

The  derivative  YP(I)  is  then  computed  as  follows: 


YP(I)  = 2*  C1*X(I)  + C2. 

A flow  chart  for  this  routine  is  shown  in  Figure  4. 


Error  Diagnostics: 
Subroutine  Required: 
Argument  List: 


Subroutine  Length: 


NONE 

QLSQ 

X = array  of  independent  variables 

Y = array  of  dependent  variables 
(displacement  or  velocity) 

YP=  array  of  derivative  data 

N = number  of  original  displace- 
ment versus  time  data  points 

NP=  number  of  points  used  to  com- 
pute each  derivative  point 

ID  = 1 — array  Y contains  displace- 
ment data  and  array  YP  will 
contain  velocity  daua 

ID  = 2 — array  Y contains  velocity 
data  and  array  YP  will  contain 
acceleration  data 


1178 
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Call  ALSQ  to  Compute 
Coefficients  for 
Point. 


Compute  I*"*1  Derivative 
Value : 


YP(I)=2*C,  XX (I)  + C2 


2.1.4  Subroutine  QLSQ  (X,  Y,  Nl,  N2 , C) 

Subroutine  QLSQ  uses  the  method  of  least  squares  to 
compute  the  quadratic  coefficients  (C^,  C2 , and  C^)  for  an 
equation  of  the  form: 

Y = C1*X2  + C2  * X + C3 


for  FN  data  points  (FN  = N2  - Nl  + 1)  from  X and  Y array  indices 
Nl  to  N2 . FN  must  be  an  odd  integer  >.  3. 

Method 

The  independent  variable  X(I)  is  transulated  by  a factor 
FF,  where 


FF  = X (NN) , 


NN 


Nl  + N2 
2 


and  XP(I)  = X(I)  - FF. 

The  quadratic  equation  in  terms  of  the  transulated  independent 
variable  is 


Y = A^  * XP2  + A2  * XP  + A3> 

The  least  square  residuals  are  a minimum  when  the  following  equations 
are  satisfied: 

A±  * Z XP4  + A2  * E XP3  + A3  * Z XP2  = Z XP2*  Y 

A1  * Z XP3  + A2  * Z XP2  + A3  * Z XP  = Z XP*Y 

A1  * Z XP2  + A2  * Z XP  + A3  * FN  = EY 

where  summations  of  XP  and  Y are  computed  for  index  I equal  Nl  to 

N2.  Determinants  are  used  to  solve  the  above  system  of  equations 
for  the  coefficients  A^  A2 , and  A3«  The  , C2 , and  C3  coefficients 
are  computed  from  A^,  A2 , and  A3  as  follows: 


I 
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r 


C2  = A2  - 2*  Al  *FF 

C3  = A3  + Ax  * FF2  - A2  * FF. 


A flow  chart  for  this  routine  is  shown  in  Figure  5 . 

Error  Diagnostics:  NONE 

Subroutines  Required:  NONE 

Argument  List:  X = array  of  independent  variables 

Y = array  of  dependent  variables 

Nl  = index  of  first  point  used 
in  fit 

N2  = index  of  last  point  used 
in  fit 

C = array  containing  quadratic 
coefficients • 


2.2  DATA  PREPARATION  FOR  INPUT  TO  HIFPD 

Preparation  of  data  for  input  to  HIFPD  consists  of 
editing  and  digitizing.  The  editing  function  provides  film 
f rame-to-time  conversion  and  film  frames  coordinates  to  plane  of 
motion  coordinates  conversion  factors.  The  digitizing  function 
provides  the  frame-by-frame  "reading"  of  the  projected  film 
frame  coordinates.  The  references,  or  "standards,"  required  to 
process  the  data  are  film  time  reference  pulses  and  surveyed 
fiducials  in  two  planes  normal  to  the  optical  axis  of  the 
camera. 

Timing  of  the  film  frames  was  accomplished  by  calculating 
the  average  film  speed  over  a span  of  approximately  150  frames. 
The  procedure  employed  to  determine  the  average  film  speed  is 
described  in  Appendix  B,  Paragrah  3.2. 
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Figure  5.  QLSQ  Flow  Chart. 
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The  first  frame  in  which  the  stroboscopic  flash  was 

observed  was  defined  as  t = 0.  The  strobe,  initiated  by  a time 

synchronizing  pulse  which  was  also  recorded  on  the  magnetic  tape 

recordings,  actually  gives  tQ  indication  within  2.0  milliseconds 

accuracy  at  the  nominal  film  speed  of  500  frames  per  second 

with  a 140°  shutter.  Since  the  flash  is  not  observed  in  film 

frame  -0001  and  is  observed  in  film  frame  0000,  it  is  apparent 

that  it  was  initiated  between  the  closing  of  the  shutter  on  film 

frame  -0001  and  the  closing  of  the  shutter  on  film  frame  0000. 

During  most  tests,  the  intensity  of  the  first  observed  flash 

would  indicate  that  it  was  initiated  between  the  closing  of  the 

shutter  on  frame  -0001  and  the  opening  of  the  shutter  on  frame 

000.  If  this  is  the  case,  the  t indication  could  be  considered 

o 

to  be  accurate  to  1.2  milliseconds,  i.e.,: 


360°  - 140°  „ , „„ 

ocTTo — x 2 msec  = 1.22  msec. 


Determination  of  conversion  constants  to  be  applied  to  the 
digitized  readings  of  the  anthropometric  points  on  the  subject 
required  that  the  following  be  known. 

a.  The  distance,  normal  to  the  plane  of  symmetry  of  the 
subject,  from  that  plane  to  the  left  edge  of  the  seat  pan. 

b.  The  distance,  normal  to  the  plane  of  symmetry  of 
the  subject  from  that  plane,  to  the  lexan  panel  above  the  side 
edge  of  the  sled  deck. 

c.  The  distances,  normal  to  the  plane  of  symmetry  of 
the  subject,  from  that  plane  to  the  anthropometric  points  to  be 
tracked. 

d.  That  the  optical  axis  of  the  primary  camera  was 
normal  to  the  plane  of  symmetry  of  the  subject. 

e.  The  distances,  between  centers,  of  the  fiducials 
mounted  on  the  edge  of  the  seat  pan  and  of  the  fiducials  mounted 
on  the  lexan  panel. 
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The  coordinates  of  the  reference  fiducials  on  the  seat 
pan  and  the  lexan  panel  were  digitized  five  times.  The  readings 
of  these  coordinates  were  then  averaged  and  the  digital  distance 
between  the  averaged  coordinates  of  each  pair  was  calculated. 
Dividing  each  of  these  digital  distances  by  the  corresponding 
measured  dimension  between  fiducials  yielded  conversion  constants, 
in  terras  of  "counts  per  foot",  in  two  planes  normal  to  the  optical 
axis.  Having  determined  these  conversion  constants,  and  having 
measured  the  distance  between  the  parallel  planes  in  which  the 
fiducials  lay,  the  distance  along  the  optical  axis  from  the 
focal  point  of  the  lens  to  each  of  these  planes  and  the  plane 
of  symmetry  could  then  be  calculated. 

Prior  to  each  test  run  the  breadth  of  the  subject  was 
measured  at  each  tracking  fiducial  location  with  an  anthropometer . 
Assuming  that  each  subject  was  symmetrical,  the  distance  from 
the  plane  of  symmetry  to  each  tracking  fiducial  was  defined  as 
one-half  the  measured  breadth  of  the  subject  of  each  fiducial 
location.  Conversion  constants  for  each  plane  parallel  to  the 
plane  of  symmetry,  thus  normal  to  the  optical  axis  in  which  a 
tracking  fiducial  lay  were  then  calculated  by  similar  triangles. 

The  actual  digitization  of  the  photometric  data  was 
accomplished  on  a Producers  Service  Corporation  model  PVR  film 
analyzer.  The  magnification  factor  of  the  projector  was 
approximately  19.6,  giving  a projected  frame  image  of  8 x 5.8- 
inches.  The  optical  encoders  were  coupled  to  the  reading  cross- 
hairs through  drive  wires  in  such  a manner  that  a displacement 
of  either  crosshair  of  one-inch  caused  the  associated  optical 
encoder  to  increase  or  decrease  the  reading  by  one  thousand 
count. 

The  operator  located  the  first  frame  in  which  the 
synchronizing  flash  was  observed  and  reset  the  frame  counter 
to  0000.  The  optical  center  of  the  film  frame  was  found  by 
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numerically  bisecting  the  vertical  and  horizontal  dimensions  of 
the  frame  image.  The  operator  then  positioned  the  crosshairs 
over  the  range  reference  fiducial  and  depressed  the  record  switch 
causing  the  frame  number  and  coordinates  of  the  fiducial  to  be 
punched  into  paper  tape  and  typed  on  the  carriage  of  the  teletype 
terminal.  He  then  proceeded  to  position  the  crosshairs  over  the 
seat  reference  fiducial.  Again,  depressing  the  record  switch 
caused  the  coordinates  to  be  recorded  on  the  listing  and  the 
paper  tape.  In  this  manner  he  would  proceed  to  each  of  the  other 
points  mentioned  in  Section  2.1,  recording  their  coordinates, 
until  all  readings  had  been  extracted  from  that  frame. 

After  advancing  the  film  to  the  next  frame,  the  operator 
would  check  the  coordinates  of  the  range  and  seat  fiducials.  If 
the  frame- to- frame  variation  of  these  coordinates  exceeded  ± 10 
counts  he  would  again  locate  the  optical  center  of  the  film 
frame  image  before  proceeding. 

This  procedure  was  repeated  for  each  film  frame  until 
the  subject  appeared  to  have  attained  a static  position  after 
the  impact. 

The  resulting  paper  tape  was  read  into  file  on  the 
CYBER  74  at  the  ASD  computer  center  from  the  teletype  terminal 
via  data  modem  on  a voice  quality  line,  and  the  file  was  edited 
and  corrected  when  necessary.  At  this  time  the  control  and 
title  data  were  added  to  the  file.  This  file  was  then  copied 
on  the  card  punch  and  printer  as  a time  saving  measure  in  case 
the  disk  file  should  be  accidentally  purged. 

At  this  point  the  program  HIFPD  could  have  been  attached 
and  executed;  however,  the  normal  procedure  was  to  obtain  the 
card  files  and  submit  them  in  the  batch  mode  on  an  overnight 
schedule.  This  permitted  the  connect  time  to  be  used  for 
read-in  and  editing  of  additional  data  files. 

A typical  standard  practice  procedure  detailing  the 
data  preparation  is  given  in  Appendix  B. 
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SECTION  3 

ANALYSIS  OF  DATA  FROM  TWO  CAMERAS 

Tracking  of  objects  moving  through  the  space  of  a three- 
dimensional  coordinate  system,  using  analysis  of  phototheodolite 
recordings,  has  been  accomplished  with  a great  degree  of  accuracy 
for  many  years.  The  recognized  requirements  to  accomplish  these 
analyses  are  quite  simple,  viz.,  cameras  at  two  or  more  tracking 
stations  with  precisely  synchronized  shutters,  time  for  frame 
signal  numbers,  azimuth  and  elevation  indicators,  and  relatively 
massive  mounted  systems  having  provision  for  leveling  and  plumbing 
the  cameras.  Additionally,  the  locations  of  the  camera  mounts 
must  be  accurately  surveyed. 

When  one  attempts  to  track  the  motion  of  several  points 
with  respect  to  a three-dimensional  coordinate  system  on  a moving 
vehicle,  all  the  simplicity  of  the  large  tracking  range  disappears, 
and  an  entirely  new  set  of  problems  confront  the  investigator. 
Camera  weight  must  be  kept  to  a minimum  and  light-weight,  rigid 
mounting  systems  must  be  devised.  Actual  shutter  synchronization 
is  not  practical  at  the  high  film  transport  speeds  necessary  to 
record  the  motion.  The  precise  locations  and  orientations  of 
the  cameras  are  subject  to  change,  not  only  from  test  to  test, 
but  even  during  a test. 

To  solve  the  above  mentioned  problems,  two  processes 
were  developed,  one  to  accurately  determine  the  location  and 
orientation  of  a camera  with  respect  to  the  coordinate  system, 
and  the  second  to  track  the  coordinates  of  identified  points 
with  respect  to  the  origin  of  the  coordinate  system. 

The  processing  programs  developed  to  accomplish  these 
solutions,  "POOCH"  and  "SLED"  respectively,  are  described  in 
the  following  discussions. 


3.1  CAMERA  LOCATION  AND  ORIENTATION 

The  value  of  a program  which  can  determine  the  location , 
orientation,  and  focal  length  of  a camera  by  examination  of 
the  film  frame  images  of  a set  of  accurately  surveyed  object 
prints  is  rather  obvious.  Regardless  of  the  care  exercised, 
if  a large  enough  amount  of  camera  data  is  taken,  sooner  or  later 
a camera  will  be  jarred  or  moved  or  a set  of  camera  data  will 
be  misplaced  or  mislabeled,  or  a camera  mounted  aboard  a moving 
vehicle  may  vibrate  on  its  mount  and  one  may  wish  to  do  a frame 
by  frame  correction  of  the  camera  position  and  orientation. 

There  are  seven  unknows  to  be  determined  in  this 
problem.  They  are  the  location  of  the  camera  focal  point  (xx,  yy, 
zz) , the  camera  focal  length  (ff)  , and  azimuth  angle  (th)  and 
the  elevation  angle  (ph)  of  the  camera  optical  axis,  and  the 
angle  (ct)  by  which  the  camera  is  tilted  from  horizontal  about 
its  optical  axis. 

Before  proceeding  with  the  main  theme,  a comment  must 
be  made  regarding  the  relationship  between  camera  focal  length 
and  the  distance  of  the  camera  from  the  scene.  There  is  a 
misconception  that  if  a camera  of  one  focal  length  is  used  to 
photograph  a scene,  a camera  of  a shorter  focal  length  will 
photograph  precisely  the  same  scene  if  taken  somewhat  closer. 

That  is  not  true  in  general.  It  is  true  only  in  the  case 
that  every  point  of  the  scene  being  photographed  lies  in  a 
plane  normal  to  the  optical  axes  of  the  two  cameras. 

Although  not  physically  correct,  one  loses  nothing 
mathematically  by  assuming  that  the  film  frame  is  interposed 
between  the  camera  focal  point  and  the  scene  and  that  film 
frame  images  are  created  by  light  rays  emanating  from  the 
object  points  and  passing  through  the  film  to  the  camera 
focal  point.  See  Figure  6. 
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Figure  6. 

first  that  a camera  of  focal  length  ff  lies  at  a dis- 
d from  a scene  composed  of  two  object  points  at  heights 
hb  above  its  optical  axis.  See  Figure  7. 


Figure  7. 


Both  object  points  lie  in  a plane  normal  to  the  camera  optical 
axis.  By  similar  triangles 


If  k > 0,  by  the  equation 

kd  d 

kff  “ ff  ' 

It  is  obvious  a camera  at  distance  kd  from  the  object  points  with 
a focal  length  of  kff  gives  exactly  the  same  image  points  at 
the  original  camera  at  distance  d with  focal  length  ff.  Thus, 
if  the  object  points  are  all  in  a plane  normal  to  the  camera 
optical  axis,  it  is  impossible  to  determine  from  the  image  and 
object  point  locations  both  the  camera  focal  length  and  its 
distance  from  the  scene. 

If  the  two  object  points  have  different  spacings  along 
the  camera  optical  axis,  the  situation  is  quite  different.  See 
Figure  8. 


Figure  8 . 


In  Figure  8 the  project  of  the  two  object  points  on  the  camera 
optical  axis  are  spaced  a distance  e apart  where  e is  known 
and  e --0.  da,  db,  ha,  and  hb  are  known.  ff  and  d are  unknown. 
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Again  by  similar  triangles 


From  (1) 


Substituting  in  (2) 


(1) 

d 

ha 

ff 

da 

(2) 

d + e 
hb 

ff 

db 

, ha 

d " da 

ff 

ha  f f 

e 

da 


hb 


ff 

db 


or 

ha  ff  + da  e 

ff 

da  hb 

" db 

or 

db  ha  ff  + da 

db  e = da  hb  ff 

(3) 

(da  hb  - db 

ha)  ff  = da  db  e 

Now  if  it  is 

assumed  that  da 

hb  - db  ha 

= 0 , we  get 

db  ha 

db  = hb  ' 

Substituting 

for  da  in  (1)  gives 

d 

ff 

hb  ff 

ha 

, db  haN 

1 hb  ' 

db  ha  ' 

or 

ff 

ha  d 

d 

db 

hb  ha 

hb  * 

But  (2)  is 

ff 

d + e 

db 

hb  ' 

Therefore , 

ff 

db 

d + e 
hb 

10  IS 

II 

e 

hb 

= 0 

e 

= 0 contrary  to  the  problem 

Therefore  in  (3)  the  coefficient  of  ff  / 0 and  we  may  solve 

for  ff: 


and  we  have  in  (4)  and  (5)  unique  solutions  for  both  ff  and  d. 

However,  two  object  points  with  different  spacings  along 
the  camera  optical  axis  are  not  sufficient  to  locate  the  camera 
becasue  the  camera  may  be  rotated  in  space  about  an  axis  between 
two  object  points  without  disturbing  any  of  the  geometrical 
relationships  in  Figure  8.  See  Figure  9. 


0) 


Figure  9. 


Figure  9 shows  the  axis  of  rotation  u>  and  the  radius  r 
of  the  circle  the  camera  focal  point  can  describe  in  space.  A 
third  object  point  lying  on  the  axis  between  the  first  two 
object  points  does  not  help  to  locate  the  camera.  To  be  of  help, 
the  third  object  point  must  lie  at  a distance  from  the  line  passing 
through  the  other  two  object  points.  The  conclusion  is  then  that 
in  order  to  completely  locate,  orient,  and  determine  the  focal 
length  of  a camera,  one  requires  an  absolute  minimum  of  three 
accurately  surveyed  object  points  forming  a nondegenerate 
triangle  whose  plane  was  not  normal  to  the  camera  optical  axis 
when  the  picture  was  taken. 
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With  the  minimum  of  three  object  points,  small  reading 
errors  can  easily  distort  the  solution.  In  practice  it  is  best 
to  have  at  least  five  to  eight  object  points  with  a good  spread 
normal  to  the  camera  optical  axis  and  significant  differences  in 
spacing  along  the  optical  axis. 

The  above  discussion  demonstrating  the  necessity  of  at 
least  three  object  points  is  in  no  way  indicative  of  the  solution 
scheme  used  in  the  camera  location  program.  The  analysis  used 
in  the  program  will  now  be  discussed. 

3.1.1  Mathematical  Models 

The  problem  is  modeled  or  indicated  in  Figure  10.  The 
principal  coordinate  system  is  a right  hand  mutually  orthogonal 
x y z coordinate  system.  Unit  vectors  e,  n,  and  h are  directed 
along  the  positive  x,  y,  and  z axis  respectively.  The  camera 
position  and  orientation  will  be  determined  with  respect  to  the 
x,  y,  z coordinate  system.  There  are  n accurately  surveyed 
object  points  at  (xi,  vi,  zi) , i = l,n.  The  camera  itself  is 
modeled  as  a film  frame  interposed  between  the  camera  focal 
point  at  (xx,  yy,zz)  and  the  scene  (consisting  of  n object  points), 
fn  is  a unit  vector  along  the  camera  optical  axis  normal  to  the 
film  frame,  p and  q are  mutually  orthogonal  film  frame  axes. 

Unit  vectors  T and  "j  are  directed  along  the  positive  p and  q 
axes  respectively.  The  p and  q axes  are  respectively  horizontal 
and  vertical  with  respect  to  the  film  frame.  Throughout  this 
report  ff  will  be  referred  to  as  the  camera  focal  length,  but 
it  is  actually  the  product  of  the  true  camera  focal  length  and 
the  magnification  of  the  film  frame  reader.  To  each  of  the 
object  points  (xi,  yi,  zi) , i=l,  n corresponds  to  a film  frame 
image  point  with  known  film  frame  coordinates  (pi,  qi) , i=  l,n. 

The  unit  vectors  fn,  T,  and  5 form  a mutually  orthogonal  right 
hand  set,  i.e.,  J = fn  x T. 

To  reiterate,  the  unknowns  in  this  problem  are  the 
location  of  the  camera  focal  point  (xx,  yy,zz),  the  focal  length 
ff,  the  azimuth  angle,  th,  and  the  elevation  angle,  ph,  of  the 
camera  optical  axis,  and  the  angle,  y , by  which  the  camera  is 
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Figure  10 


tilted  from  horizontal  about  its  optical  axis.  Of  these  seven 
unknowns,  the  last  three  are  not  explicit  in  Figure  10.  If  the 
components  of  the  unit  vector  fn  are  found  with  respect  to  the 
x y z system,  then  the  azimuth  and  elevation  angles  of  the  camera 
optical  axis  are  easily  found.  See  Figure  11. 
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Figure  11. 

fn  = cos  ph  cos  th  e + cos  ph  sin  th  n + sin  ph  h 
= fnl  e + fn2  n + fn3  h,  if  ph  and  th  are  known. 

Similarly,  if  fnl,  fn2,  fn3  are  known,  i.e., 

fn  = fnl  e + fn2  n + fn3  h 


Then 


ph  = sin  1 (fn3) 

th  = cos'1  f yj  fnl2  + fn22) 

th  = -th  if  fn2  < 0. 


Therefore,  knowledge  of  Tn  is  precisely  equivalent  to  knowledge 
of  th  and  ph  and  vice-versa.  Similarly  T determines  y but  the 
details  will  be  deferred  until  later. 

The  solution  is  done  in  two  parts.  In  the  first  part 

the  unknowns  xx,yy,zz,  and  ff  are  found.  The  angle  subtended 

at  the  camera  focal  point  by  the  ith  and  jth  image  points  is 
• • 

gy  (j  f i) . See  Figure  12. 
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Figure  12. 

Although  the  components  of  the  unit  vectors  fn,  i,  and  J are 
not  known  at  present  with  respect  to  the  xyz  system,  it  is  known 
that  they  form  a right  hand  mutually  orthogonal  set.  So  vectors 
fT  and  fj  to  the  ith  and  jth  image  points  may  be  written: 

fT  = ff  fn  + pi  X + qi  j 

f j = ff  fn  + pj  r + qj  j 


Then  the  angle  gy  subtended  by  fi  and  f j at  the  camera  focal 
point  is : 


age  point 
oordinate 
le  hy  su 
x,  yy,  zz 
ows . See 


aring  in  the  equation  for  gy  is  ff. 
tes  (pi,  qi) , i = 1,  n and  all  object 
, zi) , i = 1,  n are  known).  Similarl 
t an  arbitrary  chosen  focal  point  loc 
ith  and  jth  object  points  is  obtained 


l(xi-xx)  + (yi-yy)  + (zi-zz)  ] [(xj-xx)  +(yj-yy)  +(zj-zz) 
and  the  only  unknowns  in  the  hy* equation  are  xx,  yy,  and  zz. 


Thus  we  have 


gy  = gy  (ff)  , 


a function  of  ff,  and 

• • • • 

hy  = hy  (xx,yy, zz) , 

a function  of  focal  point  position.  If  the  correct  values  of  ff, 
x*, yy,  zz  are  inserted  in  these  two  equations,  we  should  have 

• • •• 

gy(ff)  = hy  (xx,yy,zz). 

• • ••  • • 2 

. Define  erry  = (gy  - hy)  , and 

erry  = erry  (ff ,xx,yy,zz) . 

Then  the  total  error  err  taken  over  ail  possible  pairs  of  points 
is,  for  an  arbitrary  choice  of  ff,  xx,  yy,  zz 

n-1  n 

err  = Z Z errij  (=err(ff ,xx,yy,zz) ) . 
i=l  j-1 

Then  we  should  have  at  the  correct  values  ff,  xx,  yy,  zz  err  = 0. 

Due  to  the  film  frame  reading  errors,  errors  in  the  surveyed 
object  points,  lens  distortion,  err  ^ 0 at  the  correct  value  of 
ff,  xx,  yy,  zz,  but  it  will  attain  a relative  minimum  there. 

Let  fp  and  fq  be  numbers  chosen  so  as  to  be  generous  bounds  on 
the  true  focal  length  ff.  0 < fp  <ff  < fq.  Let  (xq,yq,zq)  be 
the  x y z values  that  minimize  err  (f q , x,  y,  z) , that  is,  err 
(fq  , xq,  yz,zq)  is  the  smallest  value  of  err  that  can  be  obtained 
for  focal  length  of  fq.  Similarily,  let  err  ( fp , xp,yp , zp)  be 
the  minimum  of  err  for  focal  length  fp.  As  one  would  expect, 
the  optimized  position  of  the  focal  point  is  closer  to  the  scene 
for  the  smaller  value  of  focal  length.  See  Figure  14. 
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Figure  14. 

Let  fs  be  any  focal  length  value  such  that  f p < fs < fq.  Let 
(xs,  ys,  zs)  be  that  point  that  minimizes  err,  i.e,  err  (fs, 
xs,yx,  az)  is  a minimum  for  focal  length  fs.  Then  the  following 
proportions  furnish  excellent  estimates  of  xs,  ys,  and  zs  for 
a given  focal  length  f s . 


xs 

-_*P  - 

fs  - 

xq 

- xp 

fq  - 

fp 

ys_ 

fs  - 

_ £e_ 

yq 

- yp 

fq  - 

fp 

zs 

-_ZP 

fs  - 

zq 

- zp 

fq  - 

fp 

as 

fs  varies 

from 

fp  to  fq,  it  is  very  nearly 

true  that  the  point  (xs,  ys,  zs)  moves  a proportional  amount 
of  the  distance  from  p to  q along  the  line  segment.  Solving 
the  equation  for  (xs,  ys,  zs)  gives 
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xs  = xp  + (xq  - xp) 


fs  - ft 


fq  - fp 


ys  = yp  + (yq  - yp) 


zs  = zp  + (zq  - zp) 


fs  - ft 


fq  - fp 


fs  - fp 

fq  - fp  ' 


The  application  of  these  equations  will  be  discussed  shortly. 

A graph  of  the  minimized  value  of  err  versus  focal  length 
has  the  form  illustrated  in  Figure  15. 


Figure  15. 

The  relative  minimum  of  err  for  fs  [fp,fq]  occurs  at  ff.  Of 
course  as  fs  + ® , both  the  scene  and  the  film  frame  shrink  to 
points  with  respect  to  the  camera  focal  point,  so  err  will 
approach  the  fs  axis  asymptotically  from  above. 

When  data  is  input  to  the  program,  rough  estimates  of  ff, 
xx,  yy,  zz  are  read  in.  Also  input  are  two  generous  bounds  on 
ff,  fp,  and  fq  where  0 < fp  < ff  < fq.  The  first  step  is  to  locate 
the  points  P and  Q.  Assume  we  are  searching  for  point  P in 
Figure  14,  that  point  (xp,  yp,  zp)  that  makes  err  (fp,  xp,  yp, 
zp)  a minimum  for  focal  length  fp. 
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A cube  of  27  points  is  centered  at  the  input  point 
(xx,  yy,  zz)  . See  Figure  16.  The  cube's  edges  are  parallel  to 
the  principal  coordinate  system  axes  and  it  has  edge  length  2 dp 
where  dp  is  an  internally  defined  value.  Err  (fp,  x,  y,  z)  is 
computed  at  each  of  the  cube  points.  Then  (xx,  yy,  zz)  are  set 
equal  to  the  coordinates  of  that  point  that  resulted  in  the 
smallest  value  of  err.  Then  a new  cube  is  constructed  centered 
on  a new  point  (xx,  yy,  zz)  and  the  process  repeated.  Ultimately 
a cube  is  found  whose  center  point  gives  the  smallest  value  of 


Figure  16. 


err.  This  cube  then  contains  the  optimum  point.  dp  is  then 
reduced  to  a tenth  of  its  previous  value,  a new  smaller  cube 
centered  on  the  minimum  point  is  constructed  and  the  whole 
process  is  repeated  again.  This  iterative  process  is  terminated 
when  dp  falls  below  a certain  test  va^ue.  Then  (xp,  yp,  zp)  is 
set  equal  to  the  coordinates  of  the  last  minimum  point  deter- 
mined. This  is  the  process  used  to  determine  for  ff  = fp,  and 
ff  = fq  the  best  camera  focal  point  locations  (xp,  yp,  zp)  and 
(xq,  yq,  zq)  respectively.  The  above  process  is  capable  of 
moving  the  focal  point  many  feet  through  space  from  a grossly 
inaccurate  first  estimate  to  the  final  optimizing  point.  But 
this  iteration  is  quite  time  consuming.  For  example,  if  20  cubes 
are  required  to  locate  an  optimum  point,  it  requires  the  com- 
putation of  err  at  20  x 27  = 540  points.  And  if  there  are,  say, 

9x8 

9 object  points,  there  are  — 2 = 36  values  of  hy'  and  gy' 

that  must  be  computed,  differenced,  squared,  and  added  for  each 
computation  of  err.  For  this  reason,  the  above  process  is  not 
used  to  make  the  final  determination  of  the  focal  length.  Rather 
for  each  value  of  focal  length  fs  used  in  this  determination,  a 
very  nearly  optimum  point  for  optimizing  err  is  gotten  from  the 
three  proportions  mentioned  earlier  on  page  38. 

xs  = xp  + z ~Ep  <xc?  ~ XP) 

ys  = yp  + jq-  l -fP  (yq  - yp) 

ZS  = zp  + £q-  ~ fP  (zq  - zp). 

Then  err  is  minimized  at  fs  by  checking  points  very  close 
to  (xs,  ys,  zs)  by  a simpler,  faster  process  to  be  described 
shortly.  To  find  the  correct  value  of  ff  we  set 

dist  = (fq  - fp)  / 6.0. 

Then  fs  is  set  equal  to  fq  + i dist,  i = 0,  1,  2,  ...,  6 in 
succession  and  the  minimum  error  values  err  found  for  each  fs. 
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Say  it  occurs  at  fs  = fso.  Then  we  redefine  fq  = fso  + dist 
and  fp  = fso  - dist,  set  dist  = (fq  - fp)/  6.0  and  go  through 
the  same  process  all  over  again.  This  is  continued  until  dist 
falls  below  a present  value.  Then  ff  is  set  equal  to  the  last  fs 
computed  and  (xx,  yy,  zz)  to  the  last  optimized  value  of  (xs,  ys, 
zs)  . This  is  the  final  step  in  solving  for  these  first  four  of 
the  seven  unknowns  in  this  problem. 

In  order  to  construct  the  routine  which  is  used  to  optimize 
the  initial  values  of  (xs,  ys,  zs)  gotten  from  the  three  equations 
above,  reference  must  be  made  to  Figure  17. 


u 


Figure  17. 

a = fp  e + yp  n + gp  h 
b = xq  e + yq  n + zq  h. 

The  unit  vectors  defined  as  follows  are  the  directions  in  space 
along  which  the  iterative  procedure  searches  in  order  to  optimize 
the  initial  values  of  (xs,  ys,  zs)  to  find  the  minimum  value  of 
err  for  focal  length  fs. 
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/h  x s/ 

U = S X t . 

Defined  in  this  way,  s,  t,  and  u form  a right  hand  set  of 
mutually  othogonal  unit  vectors  with  s lying  along  the  line 
segment  between  points  p and  q.  A maximum  of  six  computations 
of  err  are  required  for  each  iteration  as  opposed  to  the  27 
required  for  the  cube  search.  The  order  in  which  the  points  are 
tried  is  indicated  in  Figure  18. 


5 u 


Figure  18. 


The  value  of  err  is  known  at  the  junction  of  s,  f:,  and  u (referred 
to  as  the  center  point).  The  points  1 through  6 are  tried  in  order. 

If  the  value  of  err  at  any  of  points  1 through  6 is  smaller  than 
that  at  the  center  point,  the  next  step  will  be  in  the  same  direction. 
For  example,  if  point  3 should  be  the  first  point  to  yield  a smaller 
value  of  err  than  that  at  the  center  point,  the  routine  will  continue 
to  step  in  the  t direction  so  long  as  successively  smaller  values 
of  err  are  obtained. 

Similarly  as  for  the  cube  search,  when  all  of  points  1 through 
6 yield  a higher  value  of  err  than  that  at  the  center  point,  the 
step  size  is  reduced  and  the  whole  process  repeated.  This  continues 
until  the  step  distance  falls  below  a certain  present  value. 

The  location  of  the  camera  focal  point  and  its  focal  length 
have  now  been  found  to  be  (xx,  yy,  zz)  and  ff  respectively.  As 
mentioned  previously,  the  determination  of  the  components  of  the  unit 
vectors  fn  and  i with  respect  to  the  principal  coordinate  system  is 
equivalent  to  finding  the  azimuth  and  elevation  of  the  optical  axis 
and  the  tilt  angle  of  the  camera  from  horizontal  about  the  optical 
axis  (refer  to  Figure  10). 


The  whole  point  of  the  following  is  to  determine  with  respect 
to  the  principal  coordinate  system  the  approximate  location  of  the 
center  of  the  film  frame,  the  origin  of  the  p,  q film  frame  coordinate 
system.  Three  image  points  and  their  respective  object  points  are 
required  to  get  a solution  and,  in  order  to  minimize  the  effects 
of  reading  and  other  errors,  they  are  chosen  in  such  a way  as  to 
form  a large  triangle  on  the  film  frame.  To  this  end,  three  sorts 
are  performed  on  the  image  points.  (pk  qk)is  that  point  farthest 
from  the  film  frame  origin.  (pi,  ql)  is  that  point  farthest  from 
(pk,  qk) . (pm,qm)  is  that  point  the  sum  of  whose  distances  from 
(pk,qk)  and  pl,ql)  is  the  greatest.  See  Figure  19. 

The  vectors  from  the  camera  focal  point  to  the  kth,  1th, 
and  mth  image  points  are  respectively: 
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skc  = 

ff 

fn 

+ 

pk 

i 

+ 

qk 

j 

sic  = 

ff 

fn 

+ 

Pi 

T~ 

1 

+ 

qi 

J 

smc  = 

ff 

fn 

+ 

pm 

-T“ 

1 

+ 

qm 

— 

D 

q 


Figure  19. 

where  these  three  vectors  have  been  written  in  terms  of  their 
components  with  respect  to  these  as -yet- unknown  unit  vectors,  fn 
i,  j of  the  camera  coordinate  system.  The  vector  magnitudes  are 


dfk  = (ff2 

+ 

pk2 

+ 

2 X/2 
q*  ) 

dfl  = (ff2 

+ 

Pi2 

+ 

2 1/2 
qi  ) 

dfm  = (ff2 

+ 

2 

pm 

+ 

2 1/2 
qm  ) 

See  Figure  20. 


(pk , qk) 


Figure  20. 


If  the  camera  was  properly  oriented,  then  skc,  sic,  sle,  and  sine 
should  point  toward  the  kth,  1th,  and  mth  object  points  respectively 
Unit  vectors  nfa,  nfb , and  nfc  pointing  at  the  kth,  1th,  and  mth 
object  points  respectively  obtained  as  follows: 


fa  = (xk  - xx)  e + 

fb  = (xl  - xx)  e + 

fc  = ( xm  - xx)  e + 

dam  = /fa/  = [ (xh  - 
dbm  = /fb/  = [(xl  - 
dcm  = /fc/  = [(xm  - 


(yk  - yy)n  + (zk  - 
(yl  - yy)n  + (zl  - 
(ym  - yy)n  + (zm  - 
xx)2  + (yk  - yy)2  + 

xx) 2 + (yl  - yy) 2 + 

2 2 
xx)  + (ym  - yy)  + 


zz)  h 
zz)  h 
zz)  h 

(zk  - zz) 
(zl  - zz) 
(zm  - zz ) 


, 1/2 

? 1/2 
Z] 

y 1/2 
2] 


nfa  = fa/dam 
nfb  = fb/dbm 
nfc  = fe/dem. 


Were  it  not  for  the  various  errors,  the  above  equivalence  would 
be  exact  and  the  unknown  unit  vectors  fn,  T,  and  J could  be  solved 
for  directly.  The  solution  could  be  done  anyway,  but  the  resulting 
solution  vectors  fn,  i,  and  j would  in  general  be  neither  unit 
vectors  nor  would  they  be  mutually  perpendicular.  Of  course,  the 
resulting  system  could  be  made  orthonormal  algebraically,  but 
it  was  not  known  at  the  time  how  large  errors  would  be  nor  how  much 

distortion  would  be  produced  by  the  errors.  It  was  decided  to  take 

advantage  of  the  approximate  vector  equivalence  above  in  the  following 
manner  in  order  to  determine  7n,  T,  and  J. 

Let 

w = skc  - smc  = (pk  - pm)  T + (qk  - qm) J 

u = sic  - smc  = (pi  - pm)  T + (ql  - qm) J. 

A vector  c from  the  image  point  in  to  the  center  of  the  film  frame 


c = - pm  1 - qm  j 


= c a i + cb  j (defining  ca  and  cb) . 


See  Figure  22, 


Making  w and  u unit  vectors  by  setting 


w = and  u 

/w/ 


c in  terms  of  w and  u 


, w4  may  express 


c = a w + g u where  a and  £5  must  be  found. 


a wa  + 

8 

ua 

= ca 

a wb  + 

6 

ub 

= cb 

and 


Therefore,  in  terms  of  the  x y z coordinate  system,  it  is  approxi 
mately  true  that 


c 


sk  - sm 
/sic  - sm/ 


si  - sm 
/sT  - sm/ 


So  a vector  from  the  focal  point  to  the  center  of  the  film  frame 
is  approximately: 


cff  = sm  + c 


and 


fn 


cff 

/cff/ 


= fn  1 e + fn  2 n + fn  3 h. 


Figure  23. 


Then  the  azimuth  th  and  the  elevation  ph  of  the  optical  axis  are 
approximately 


I 


th  = cos 


If  (fn  2 < 0.0)  , th  = -th 


fn  1 

[(fn  l)2  + (fn  a) 


ph  = sin  (fn  3) 


(-it  < th  < ir  ) 


< i < ?h  < i > 


In  order  to  find  the  approximate  direction  of  the  film  frame  vector 
i with  respect  to  the  x y z coordinate  system  two  new  unit  vectors 
i'  and  J'  are  constructed.  These  lie  in  the  film  frame  as  do  i and 
J and  they  also  form  with  fn  a right-handed  trio  of  mutually 
orthogonal  vectors.  See  Figure  24. 


Figure  24. 

In  order  to  have  a unit  vector  IT  which  will  always  be  well 
defined  regardless  of  the  orientation  of  the  camera  optical  axis 
fn,  i'  is  derived  from  that  component  of  the  vector  joining  object 
points  k and  1 that  lies  parallel  to  the  film  frame,  that  is, 
normal  to  fn.  It  will  be  recalled  that  image  points  k and  1 had 
the  greatest  separation  of  any  point  pair  on  the  film  frame,  p is 
the  angle  to  be  solved  for  once  i'  and  j1  have  been  found. 


Figure  25  and  the  following  discussion  give  the  derivation 
of  i'  and  j"1*. 


pip  = (xk  - xl)  e + (yk  - yl)  n + (zk  - zl)  h 
= pipa  e + pipb  n - pipe  h. 

To  subtract  from  pip  a vector  parallel  to  fiT  leaving  a difference 
vector  pi  normal  to  fn  requires  the  solution  of  the  following 
equation  for  the  scalar  a. 

(pip  - a fn  ) . fn  = o 
pip  . fn  - a fn  . fn  = 0 
or,  since  /fn/  = (fn.fn)  1/2  = 1, 

a = pip  . fn. 

So  pi  = pipj-  pip.  fF  fn  is  parallel  to  the  film  frame,  i.e., 
normal  to  fn. 


P = -Ei— 

/pi/ 

V = fn  . IT. 

To  find  p the  kth  image  and  object  points  are  used  as  follows. 
See  Figure  26. 


r = pk  i + qk  j 
T = cos  p V + sin  pj1- 

7 = cos  (p  + ) IT  + sin  (p  + -j-  ) V 

= ~ sin  pIT  + cos  p j ' 
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I 


1 

fa  is  a vector  from  the  focal  point  to  the  inpage  point  (pk,  qk) . 

fl  = (ff2  + pk2  + qk2)1/2ilj^-xx)  i + (yk7ZyU..+,  (zk.-zz)  h] 

t (xk-xx) 2 + (yk-yy)2  + (zk-zz)2]1/2 

fa  = faa  e + fab  n + fac  h 

ff  = ff  fn  - ff  fnle  + ff  fn2n  + ff  fn3h 

= ffa  e + ffb  n + ffc  h 

r = fa  - ff. 

So  pk  i + qkj  = fa  - ff . 

Substituting  for  T and  j their  expressions  in  IT  and  p" 

gives 

pk  cos  pIT  + pk  sin  pjr  - qk  sin  pIT  + qk  cos  pp"  = fa  - ff 

fa  - ff  can  then  be  written  as 

fa  - ff  = (fa  - ff)  . r7"  IT  + (fa  - ff)  . p"  J1” 

Equating  components  of  i'  and  j'  gives  the  following  scalar 
equations  in  which  p is  the  only  unknown. 

pk  cos  p - qk  sin  p = (fa  - ff:)  . IT  = ela  (defining  ela) 

pk  sin  p + qk  cos  p = (Ta  - FF)  . p”  = elb  (defining  elb) . 

Solving  the  above  equations  gives  p. 

Then  T - cos  p IT  + sin  p j ' where  the  components  of  IT  and 
j 7 are  known  with  respect  to  the  x y z coordinate  system.  So  we 
get  i = gia  e + gib  n + gic  h where  the  gia,  gib,  gic  are 
obtained  by  substituting  for  IT  and  J1"  their  expressions  in  e,  n, 
and  h.  Then  the  azimuth  thi  and  elevation  phi  of  T are  found 
as  follows: 
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(gia  + gib  ) ' J 
If  gib  is  less  than  0.0,  set  thi  = - thi. 

phi  - sin  1 (gic) . 

The  unit  vectors  fn  and  i obtained  above  complete  the 
orientation  of  the  camera,  but  they  are  only  approximate.  They 
depend  on  the  values  of  th,  ph,  and  p.  The  values  of  th,  ph,  and 
p found  above  should  be  very  close  to  the  correct  values.  To  find 
that  combination  of  th,  ph,  and  p which  best  aligns  all  the  image 
points  we  construct  an  error  sum  as  follows.  Due  to  the  various 
errors,  the  point  at  which  a ray  from  the  focal  point  to  an  object 
point  pierces  the  film  frame  is  not  identical  to  the  image  point 
of  that  object  point.  This  is  illustrated  in  Figure  27  for  the  ith 
image  and  object  points. 


(xi ,yi , zi) 


Figure  27. 


6 


Let  (pa,  qa)  denote  the  film  frame  coordinates  of  the  ith  pierce 
point  and  erri  be  the  distance  between  (pa  and  qa)  and  the  ith 
image  point  at  (pi,  qi) . 


Then  erri  = [(pi  - pa) 


(qi  - qa)  ] 


err  = 


n 

E erri. 
i=l 


Any  change  in  any  of  the  three  angles  th,  ph,  or  p gives  a 
new  set  of  unit  vectors  fn,  i,  j and  causes  the  coordinates  of  the 
ith  pierce  point  (pa,  qa)  to  change.  That  combination  of  th,  ph 
and  rho  ( p ) which  results  in  the  minimum  value  of  err  is  con- 
sidered to  be  the  best  fit.  Of  course,  other  definitions  of  err 
are  possible,  the  most  obvious  being  to  define  err  as  the  sum  of 
the  squares  of  the  miss  distances  rather  than  as  their  sum,  but  it 
was  felt  that  this  would  allow  one  badly  read  image  point  or  badly 
surveyed  object  point  to  have  an  unduly  high  distorting  influence 
on  the  solution. 

Before  discussing  the  search  technique  whereby  th,  ph,  and 
rho  are  varied  in  an  orderly  manner  to  find  the  minimum  value  of 
err,  the  math  steps  leading  from  each  distinct  triple  of  th,  ph, 
and  rho  to  the  ith  error  contribution  will  be  summarized. 

fn  = cos  th  cos  ph  e + sin  th  cos  ph  n + sin  ph  h 

= fnl  e + fn2  n + fn3  h 

pip  = (xk  - xl)  e + (yk  - yl)  n + (zk  - zl)  h 

= pipa  e + pipb  n + pipe  h. 


(Refer  to  Figure  25  and  its  discussion.) 


= pip.fn  = pipa  fnl  + pipb  fn2  + pipe  fn3 
= pip  - a fn 

= (pipa  - a fnl)e  + (pipb  - a fn2)  n + (pipe 
= pia  e + pib  n + pic  h 


- a fn3)h 


_ 2 2 2 1/2 
pid  = /pi/  = (pia  + pib  + pic  ) 

IT=  = Ei|i  + Ei§n+2iCh 

/piy  Pld  Pid  Pid 

= pia  e + pib  n + pic  F (Redefining  pia,  pib,  pic  to  be 
components  of  the  unit  vector  I"1". ) 

j '=  fn  x P" 

= (fn  2 pic  - fn  3 pib)  e + (fn3  pia  - fnl  pic)n 
+ (fnl  pib  - fn2  pia)h 
= pja  e + pjb  n + pjc  h. 

From  the  discussion  following  Figure  26, 

T = cos  p i’  + sin  p J7". 

Inserting  the  equation  IT  and  and  grouping  e,  n,  and  h 
i = (pia  cos  p + pja  sin  p ) e + (pib  cos  p + pjb  sin  p)n 

+ (pic  cos  p + pjc  sin  pTh 

= tia  e + tib  n + tic  h. 

Similarly,  J = -sin  p IT  + cos  pj7^ 

= (-pia,  sin  p + cos  p)e  + (-pib  sin  p + pjb  cosp)n 

+ (-pic  sin  p + pjc  cos  p)h 
= tja  e + tjb  n + tjc  h 


, 


Referring  now  to  Figure  27, 

ti  is  the  vector  locating  the  pierce  point  (pa,  qa)  due 
to  the  ith  object  point,  ti  is  a scalar  multiple  of  IT  ti  = w si 
where  w is  unknown. 

The  projection  of  ti  on  fn  must  equal  the  focal  length  ff. 

tX.fn  = ff 


w si.fn  = ff. 


si . fn 


In  component  form 


si  — (xi  - xx)  e + (yi  - yy)  n + (zi  - zz)h 


(xi  - xx)  fnl  + (yi  - yy)  fn2  + (zi  - zz)  fn3 

A vector  r from  the  center  of  the  film  frame  to  the  pierce 
point  is  gotten  from  the  equation 


r * tT  - ff  fn. 


See  Figure  28. 
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Figure  29. 


Then 


erri  = ,'ri  - r/ 


2 2 
= [ (pi  - pa)  + (qi  - qa)  ] 


Then 


m 

err  = Z erri. 
i=l 

The  scheme  whereby  th,  ph,  and  rho  are  varied  to  minimize 
err  is  indicated  in  Figure  29. 

The  points  numbered  1 through  6 indicate  the  order  in 
which  th,  ph,  and  rho  are  varied.  If  point  1 (indicating  a positive 
increment  added  to  the  present  value  of  th)  gives  a lower  err 
than  exists  at  the  origin,  then  successive  increments  are  added 
to  th  so  long  as  the  value  of  err  continues  to  fall.  Points  1, 

3,  and  5 correspond  to  increases  in  th,  ph,  and  rho  respectively 
and  2,  4,  and  6 to  decreases.  Although  the  axes  in  Figure  24 
serve  merely  to  indicate  the  order  of  variation,  the  procedure 
is  precisely  the  same  as  that  used  to  minimize  err  (fr,  xx,  ys,zs) . 
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See  Figure  18. 


When  err  (th,  ph,  rho)  has  been  minimized,  the  final 
values  of  fn,  T,  and  J are  computed.  Also  computed  are  the  azimuth 
and  elevation  of  the  vector  T,  thi  ,and  phi  respectively. 


z 


At  this  point  the  camera  has  been  located,  oriented,  and 
its  focal  length  determined.  The  only  one  of  the  seven  unknowns 
not  yet  solved  for  is  the  tilt  angle,  y,  of  the  camera  from  the 
horizontal  about  its  optical  axis.  Knowing  the  azimuth  thi  and 
the  elevation  phi  of  the  unit  vector  i makes  (gam)  redundant,  but 
it  is  among  the  output  of  the  program,  so  its  derivation  will  be 
given . 

See  Figure  30. 

ih  is  a unit  vector  in  the  film  frame  which  is  parallel 
to  the  x y plane  (i.e.,  horizontal)  of  the  main  coordinate  system. 

Ih  = 

/h  x fn/ 

By  definition  of  the  cross  product  Ih  is  normal  to  both 
fn  and  h.  Therefore,  it  both  lies  in  the  film  frame  and  is  parallel 
to  the  xy  plane.  It  is  obviously  a unit  vector. 


1 


1 
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Y = cos  ^ (ih  • h) 

Y is  to  be  recorded  as  negative  when  T lies  below  Ih  and 
positive  otherwise.  The  following  test  determines  the  sign. 

tell  = fn  . IF  x I. 

If  (tell  < 0.0)  Y = “Y  • 

(for  the  orientations  indicated  in  Figure  25  one  could 
simply  give  y the  same  sign  as  h component  of  T,  but  the  situation 
is  a little  more  complicated  when  dealing  with  a nearly  vertical 
optical  axis  as  is  the  case  with  an  overhead  camera.) 

The  solution  is  complete.  To  give  a measure  of  how  good 
the  whole  situation  is,  the  following  check  was  devised.  See 
Figure  31. 


Figure  31. 


' 


A ray  is  extended  from  the  camera  focal  point  through  each  image 
point  (pi,  qi)  past  its  corresponding  object  point  and  the  distance 
vector  dis  of  the  ray's  nearest  approach  to  the  object  point  is 
computed.  If  all  distances  dis  so  computed  are  of  small  magnitude, 
the  fit  is  obviously  good. 

dh  = ff  fn  + pi  i + qi  j 
= dha  e + dhb  n + dhc  h 

where  the  final  form  is  obtained  by  expressing  fn,  T,  and  J in 
terms  of  their  components  in  the  principal  coordinate  system  and 
then  grouping  on  the  principal  coordinate  system  unit  vectors  e,  n, 
and  h. 


sh  = (xi  - xx)  e + (yi  - yy)  n + (zi  - zz)h 

= sha  e + shb  n + she  h 

dh  = — (a  unit  vector  in  the  director  of  the 

/dh/ 

ith  image  point  (pi,qi).) 

There  is  a scalar  coefficient,  acof,  such  that 

acof  dh  = sh  + dis . 

Dotting  through  by  dh 

acof  dh  . dh  = Ih  . dh  + dis . dh 

= Ih  . dh  (since  dis  is  perpendicular  to  dh) 

acof  = sh  . dh  since  Si  . dh  = 1. 

So,  dis  = acof  dh  - sh  and  the  miss  distance  is  dis  where 

dis  = /dis/. 
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These  miss  distances  have  characteristically  fallen  into 
a range  of  0.05-inches  to  0.22-inches  in  the  solutions  run  thus 
far.  Good  agreement  in  this  check  indicates  not  only  that  the 
camera  data  found  is  correct,  but  also  that  the  survey  of  the 
object  points  and  the  reading  of  the  image  points  were  probably 
both  quite  accurate.  In  closing  the  discussion  of  the  mathematics 
used,  it  should  be  mentioned  that  one  or  two  grossly  inaccurate 
image  or  object  point  locations  can  completely  swamp  out  the 
relative  minimum  indicated  in  Figure  15. 

3.1.2  Program  POOCH 

A brief  discussion  of  the  overall  program  operation  will 
be  given  before  going  into  input  and  output  in  detail.  A reduced- 
size-copy  of  the  program  is  included  in  Appendix  E.  Figure  32 
contains  a block  diagram  of  the  program. 

The  main  program  POOCH  reads  in  and  scales  the  data. 
Subroutine  FIND  together  with  its  auxiliary  subroutines,  HORSY, 
REFIN,  and  ANGLE,  computes  the  correct  values  of  focal  point 
position  (xx,  yy,  zz)  and  focal  length  ff.  FIND  calls  HORSY 
twice  to  locate  points  P and  Q (shown  in  Figures  14  and  17) . 

HORSY  performs  the  27  point  cube  search  displayed  in  Figure  16. 
Subroutine  ANGLE  computes  the  value  of  err  for  any  focal  length 
fs  and  focal  point  location  (sx,  ys , zs) . After  points  P and  Q 
(see  Figures  14  and  17)  have  been  found,  FIND  uses  subroutine 
REFIN  to  find  the  best  values  of  focal  length  ff  and  focal  point 
location  (xx,  yy,zz).  REFIN  performs  the  six  point  search 
illustrated  in  Figure  29 . Both  HORSY  and  REFIN  use  ANGLE  to 
compute  err  (fx,  xs , ys,  zs) . 

Program  POOCH  then  calls  subroutine  SERCH  which  provides 
approximate  values  of  the  camera  orientation  variables  th,  ph, 
and  rho.  Finally,  subroutine  DANG  is  called  to  optimize  the 
approximate  camera  orientation  produced  by  SERCH.  DANG’s  auxiliary 
subroutine  WRANG  computes  an  error  value  which  is  the  sum  of  the 
distances  of  all  pierce  points  from  their  corresponding  image 
points . 


Figure  32.  Block  Diagram  of  Program  POOCH. 

See  Figures  27  and  28  and  the  associated  discussion. 


The  following  list  contains  the  variables  which  must  be 
input  to  the  program: 


XX 

YY 

ZZ 

FF 


The  x,  y,  and  z coordinates  of  the  estimated 
camera  position. 

Input  in  feet. 


The  nominal  focal  length  of  the  camera  in  inches 
times  the  magnification  of  the  film  frame  reader 
used. 
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FP,  FP  and  FQ  are  lower  and  upper  bounds  respectively 

FQ  for  the  final  FQ  value  of  FF.  Both  FP  and  FQ  are 

in  inches . 

XSR  A factor  intended  to  account  for  distortion  when 
present.  0.0  is  used  in  all  our  runs.  It  is 
dimensionless . 

N The  number  of  image  points  and  associated  object 

points  to  be  read  in.  Usually  5 to  8 are  used. 

X (20 ) 

The  x,  y,  z coordinates  of  up  to  N = 20  object 
2 (20)  ' Points  <in  feet). 

P(20),  The  film  frame  coordinates  of  up  to  N = 20  image 
Q(20)  points.  Input  in  Q(20)  "counts,"  1000  counts/ 
inch . 

The  input  format  is  as  follows: 

READ  (5,2)  XX,  YY,  ZZ,  FF,  FP,  FQ,  XSR,  N 
2 FORMAT  ( 7F10 .0,15) 

Then,  for  I = 1,  N 

READ  ( 5 , 4 ) P (I),  Q (I) , X(I) , Y(I),Z(I) 

4 F0RMAT  (5F10.0) 

Where  each  image  point  P(I),  Q(I)  is  on  the  same 
card  as  its  associated  object  point  X(I)  , Y(I), 
Z(I)  • 

Multiple  cases  may  be  read  into  the  program.  Figure  33 
shows  a sample  data  deck  for  a case  run. 

There  is  at  present  a great  deal  of  working  output  taken 
from  the  program  but  most  of  it  was  used  for  debug  purposes  and 
is  of  little  consequence  now  and  can  be  eliminated.  Rather 
than  waste  many  pages  describing  this  debug  output,  that  part 
of  the  output  giving  the  camera  location  and  orientation  and 
the  various  checks  on  the  goodness  of  the  solution  will  be 
discussed. 
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Figure  33.  Sample  Data  Input  to  Program  POOCH. 

Figure  34  shows  a reduced-size  section  of  the  tab  con- 
taining the  important  output. 
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Figure  34.  Sample  Data  Output  from  Program  POOCH. 
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The  first  three  lines  of  the  tab  section  contain  the  check 
values  for  the  first  image  point.  The  first  pair  of  values  after 
the  image  point  number  are  the  p and  q coordinates  of  the  first 
film  frame  image  point  in  feet.  The  second  pair  of  values  are 
the  p and  q coordinates  of  the  pierce  point  due  to  object  point  1, 
that  is,  the  p and  q coordinates  of  the  point  of  the  film  that 
would  be  pierced  by  a ray  from  the  camera  focal  point  to  the  first 
(refer  to  Figure  31)  object  point.  The  last  value  is  the  distance 
in  feet  between  the  image  point  and  the  pierce  point.  It  must, 
of  course,  be  remembered  as  explained  previously,  that  the  magnifi- 
cation of  the  film  frame  reader  enlarges  the  actual  film  frame  size 
by  about  a factor  of  20.  However,  these  distances  may  be  usefully 
compared  both  with  one  another  in  a single  run,  and  also  between 
runs  to  see  if  reasonable  results  are  being  obtained.  To  determine 
how  far  apart  the  image  point  and  its  associated  pierce  point 
are  on  the  real  film  frame,  the  distance  given  could  be  divided 
by  about  20.4,  the  approximate  magnification  of  the  reader  used. 

The  last  line  of  output  for  the  first  image  and  object 
points  gives  the  object  point  number,  then  the  x,  y,  and  z 
coordinates  in  feet  of  the  object  point.  The  fifth,  sixth,  and 
seventh  values  are  the  x,  y,  and  z components , respectively , of  the 
miss  vector  dis  described  in  Figure  31  and  the  text  pertaining 
thereto.  Its  components  are  given  in  feet.  The  last  value  is 
the  magnitude  of  dis  in  feet.  This  is  the  most  important  value 
for  judging  how  good  the  solution  is.  The  numbers  written  into 
the  right  of  these  values  are  the  magnitude  of  dis  in  inches , 
and  as  can  be  seen,  are  quite  reasonable.  This  assures  that  not 
only  did  the  program  provide  a good  solution,  but  also  that  the 
object  points  were  correctly  surveyed  and  their  film  frame  images 
were  correctly  read.  In  short,  when  the  magnitudes  of  dis  are 
small  for  each  object  point,  it  provides  a comprehensive  check 
on  the  entire  data  reduction  scheme. 


74 


In  the  run  illustrated,  there  were  seven  points  available 
on  which  to  base  the  solution.  Immediately  after  the  listing  of 
the  above  data  for  the  last  of  these  points,  the  camera  solution 
data  is  output.  In  the  line  beginning  X,  Y,  Z,  F ...  the  first 
three  numerical  values  are  the  x,  y,  and  z coordinates  of  the 
location  of  the  camera  focal  point  (in  feet) . The  fourth  value 
is  the  fictitious  program  focal  length,  i.e.,  it  is  the  product 
of  the  true  camera  focal  length  and  the  film  reader  magnification 
and  is  given  in  feet.  The  fifth  and  sixth  values  are  the  azimuth 
and  elevation  angles  of  the  camera  optical  axis,  both  given  in 
radians.  The  last  value  in  this  line  of  output  is  the  angle  that 
has  been  referred  to  us  as  rho.  It  is  of  no  consequence  after  the 
camera  solution  has  been  generated. 

The  line  beginning  DAMG, . . gives  the  azimuth  and  elevation 
angles  of  the  film  frame  T axis  and  the  tilt  angle  of  the  camera 
from  horizontal  about  the  optical  axis  in  that  order.  All  angles 
are  in  radians . 

It  yet  remains  to  explain  a minor  point  regarding  coordinate 


systems.  The  experimental  coordinate  system  is  a left  hand  system. 
See  Figure  35. 


The  input  coordinates  of  all  object  points  (xi,  yi,  zi) 
and  the  estaimated  initial  camera  position  (xx,  yy,  zz)  are  input 
with  respect  to  this  coordinate  system.  The  program  changes  the 
sign  of  the  y coordinate  of  each  of  these  points  so  that  they  are 
properly  represented  for  the  right  hand  coordinate  system  used 
in  the  program.  See  Figure  36 . 


z 


Figure  36. 

This  change  was  made  so  that  it  would  not  be  necessary 
to  define  the  components  of  every  vector  cross  product  used  for 
a left  hand  coordinate  system  when  it  is  second  nature  to  most 
to  automatically  write  them  for  a right  handed  system. 

Moreover  all  the  camera  location  and  orientation  data 
is  output  with  respect  to  the  program  right  hand  coordinate 
system.  This  is  entirely  appropriate  since  the  camera  solution 
is  used  as  input  to  yet  another  program  which  also  contains  a 
right  hand  principal  coordinate  system.  Similarly,  the  coordinates 
of  the  object  points  are  also  output  with  respect  to  the  program 
right  hand  system  rather  than  with  respect  to  the  original 
experimental  left  hand  system. 

The  film  frame  reader  has  a coordinate  system  as 
indicated  in  Figure  37  with  the  horizontal  axis  positive  to  the 
right . 
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Figure  37 . 

In  order  that  the  unit  vectors  fn,  i,  j (refer  to  Figure 
12)  in  that  order  form  a right  hand  system,  the  sign  of  the 
horizontal  coordinate  of  each  image  point  is  changed.  Moreover 
the  film  frame  coordinates  are  output  with  the  signs  of  their 
horizontal  coordinates  still  reversed. 

There  is  nothing  very  involved  about  these  sign  changes, 
but  one  should  be  aware  of  what  has  been  done  if  he  attempts  to 
use  the  program.  These  sign  changes  are  clearer  from  the  tab  of 
program  POOCH  than  from  all  the  explanation  above,  and  could  be 
easily  changed  to  suit  if  one  were  going  to  use  the  program  for 
an  experiment  in  which  the  experimental  coordinate  system  itself 
were  right  handed. 


3.2  TWO-RAY  SOLUTIONS  OF  COORDINATES  OF  POINTS 

The  Camera  Location  Program  has  direct  application  in 
Program  Sled  which  performs  the  reduction  of  data  taken  during 
seat  harness  tests.  A dummy  and  an  animal,  or  two  dummies, 
are  strapped  into  seats  mounted  aboard  an  acceleration  sled. 
Also  rigidly  mounted  aboard  the  sled  are  two  high-speed  cameras 
which  photograph  the  reactions  of  the  test  subjects  as  the  sled 
accelerates  down  its  track  in  the  +x  direction.  See  Figure  38. 
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Figure  38. 


The  experimental  coordinate  system  is  a left-hand  x y z 
system  affixed  to  the  moving  sled.  All  subject  motion  recorded 
is  understood  to  be  relative  to  the  sled.  In  cases  run  so  far, 
two  targets  (plastic  disks  marked  so  as  to  be  easily  discerned 
on  the  film  frame  reader)  are  attached  to  each  test  subject. 

For  a dummy,  the  forehead  and  chin  have  been  targeted.  When  the 
subject  was  an  animal, head  accelerometer  pack  and  tip  of  the  snout 
were  chosen  as  target  points.  The  whole  purpose  of  program  Sled 
is  to  provide  displacement  versus  time  graphs  of  each  of  the 
targeted  points  on  each  of  the  test  subjects. 

The  camera  location  program  is  used  first  to  accurately 
locate  and  orient  each  of  the  two  cameras.  The  coordinates  of  a 
number  of  points  rigidly  affixed  to  the  sled  are  accurately 
measured.  Their  film  frame  image  coordinates  are  determined  for 
each  camera.  The  camera  location  and  orientation  data  obtained 
from  the  camera  location  program  is  then  input  to  the  Sled  program. 
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The  principal  coordinate  system  in  the  Sled  program  is  a 
right  hand  x y z system,  not  the  experimental  left-hand  system 
indicated  in  Figure  38.  This  change  is  accomplished  simply  by 
changing  the  sign  of  the  y coordinate  of  every  point  after  it 
has  been  input  to  Sled.  Also,  so  that  the  camera  system  axes  will 
constitute  a right-hand  coordinate  system,  the  horizontal  film 
frame  coordinate  is  reversed  in  sign.  Unit  vectors  e,  n,  and  h 
are  directed  along  the  x,  y,  z axes  respectively.  See  Figure  39. 
The  x,  y,  and  z axes  are  mutually  perpendicular. 

z 


Figure  39. 

The  two  cameras  are  designated  camera  a and  camera  b 
and  their  focal  points  are  located  by  vectors  ra  and  rb 
respectively.  See  Figure  40. 
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Figure  40. 


ra  = raa  e + rab  n + rac  h 

rb  = rba  e + rbb  n + rbc  h 

The  cameras  are  modeled  in  the  same  manner  as  in  the  camera 
location  program,  that  is,  as  a film  frame  interposed  between  the 
scene  and  the  camera  focal  point.  The  camera  a coordinate  system 
is  illustrated  in  Figure  41. 


Figure  41. 


fa  lies  along  the  camera  optical  axis  and  is  perpendicular  to 
the  horizontal  and  vertical  film  frame  unit  vectors  ai  and  aj. 

The  magnitude  of  fa  is  fa,  the  camera  a focal  length.  (As  in 
the  camera  location  program,  fa  is  actually  the  product  of  the 
true  camera  focal  length  and  the  film  reader  magnification.  The 
analogous  vectors  for  camera  b are  bT,  bj,  fb. 

The  data  that  must  be  supplied  to  completely  define 
camera  a consists  of  the  three  components  of  its  position  vector, 
raa,  rab,  rac,  the  azimuth,  elevation,  and  magnitude  of  fa  which 
are  tha,  pha,  fa,  the  azimuth  and  elevation  of  aT  which  are  thia, 
ohia.  For  camera  b the  analogous  quantities  are  in  the  same 
order:  rba,  rbb,  rbc,  fb,  thb,  phb,  thib,  phib. 

We  will  now  discuss  the  principal  parts  of  the  analysis 
in  the  various  subroutines.  The  first  topic  will  be  the  derivation 
of  the  a camera  ray  solution.  See  Figure  42. 


and  b respectively.  sa  and  sb  are  unit  vectors  along  the  ray 
direction  from  the  camera  focal  point  toward  the  film  frame 
image  of  one  of  the  moving  points  on  a dummy,  a point  whose 
position  is  to  be  determined  with  respect  to  the  x,  y,  z coordinate 
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system.  Due  to  film  frame  reading  errors,  camera  location  and 
orientation  errors,  and  also  to  the  interpolation  used  to  find 
sa  and  sE  at  the  same  instant  in  time,  the  rays  will  not  in 
general  intersect,  but  will  miss  each  other  by  a small  distance 
denoted  by  the  vector  r . r is  necessarily  perpendicular  to  both 
rays. 

The  position  of  the  moving  point  on  the  dummy  is  assumed 
to  lie  at  the  midpoint  of  r . ra,  rb,  la,  and  sb  are  known 
completely.  Although  the  vector  algebra  solution  below  does  not 
explicitly  use  the  names  of  the  components  of  these  vectors,  we 
will  list  them  for  completeness  sake: 

ra  = raa  e + rab  n + rac  h 

rb  = rba  e + rbb  n + rbc  h 

sa  = saa  e + sab  n + sac  h 

sb  = sba  e + sbb  n + sbc  h . 

The  components  for  ra  and  rb  are  input  to  program  SLED. 
The  method  whereby  the  components  of  sa  and  lb  are  obtained  is 
taken  up  immediately  after  the  following  vector  algebra  two  ray 
solution.  Solution: 


— sa  x sb 
rn  = ■ 

/sa  x sb/ 


is  a unit  vector  collinear  with  r,  that  is,  r is  perpendicular  to 
both  sa  and  sb. 

Let  r = d rn  where  d is  the  as-yet-unknown-scalar 
multiple  of  rn  which  yields  r.  (d  will  be  negative  if  the  sb 
ray  passes  over  the  sa  ray  in  a particular  case.) 
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(1)  ra  + coa  sa  = rb  + cob  sb  + r where  coa  and  cob 
are  scale  factors  to  be  solved  for. 

Substituting  for  r gives 

(2)  ra  + coa  sa  = rb  + cob  sE  + d rn. 

Dotting  through  by  rn  gives 

ra  . rn  + coa  sa  . rn  = rb  . rn  + cob  ib.rn+drn.rn 
ra  . rn  - rb  . rn  = d. 

(3)  or  d = (ra  - rb)  . rn 
Rewriting  (2)  in  an  altered  form  gives: 

(4)  coa  sa  - cob  sb  = rb  - ra  + d rn. 

Dotting  through  (4)  by  sa  gives: 
coa  - cob  sb  . sa  = (rb  - ra)  . sa 
or  (5)  coa  = (rb  - ra)  . sa  + cob  sa  . sb 

Dotting  through  (4)  by  sb  gives: 

(6)  coa  sa  . sb  - cob  = (rb  - ra)  . IE. 

Substituting  (5)  for  coa  in  (6)  gives 

[(rb  - rn)  . sa  + cob  sa  . sb  sa].  sb  - cob  = (rb  - ra)  . Ib. 


Solving  for  cob  gives: 


(7)  cob  - (rk  ~ ra)  » sa  sa  ■ sb  - (rb  - ra)  . sb 
1 - (sa  . sb)2 


Substituting  the  value  of  cob  into  (5)  gives  coa. 

A vector  sol  locating  the  solution  point  (midpoint  of  r) 
can  be  computed  by  either  of  two  equivalent  vector  sums. 


or 


sol  = rb  + cob  sb  + 0.5  r (see  Figure  9) 


sol  = ra  + coa  sa  = 0.5  r. 


The  ray  miss  distance  dis  is  gotten  from  dis  = /d/.  d in  the 
equation  below 

r = d rn 

will  be  negative  if  rn  should  have  a direction  or  sense  opposite 
to  r.  Its  sign  depends  on  the  order  in  which  the  cameras  were 
read  in  and  which  ray  is  uppermost. 

The  above  solution  is  performed  in  subroutine  SOLVE. 

In  order  that  the  unit  vectors  sa  and  ib  have  their  components 
computed  with  respect  to  the  main  coordinate  system,  the  com- 
ponents of  each  camera  coordinate  system  vector  must  be  known 
with  respect  to  the  main  coordinate  system.  Let  (pa,  qa)  be 
the  film  frame  coordinates  of  camera  a,  fa,  its  focal  length, 
tha  and  pha  the  azimuth  and  elevation  of  its  optical  axis,  thia 
and  phia  the  azimuth  and  elevation  of  the  camera  a sa  film 
frame  unit  vector.  See  Figure  43. 

Now  ra  = raa  e + rab  n + rac  h is  known,  its  components 
being  read  in. 

fa  = fa  cos  tha  cos  pha  e + fa  sin  tha  cos  pha  n 
+ fa  sin  pha  h = faa  e + fab  n + fac  h 

ai  = cos  thia  cos  phia  e + sin  thia  cos  phia  n 
+ sin  phia  FT  = aia  e + aib  n + aic  h 
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sa 


ssa 


ssa 

ssa 


where 


/ssa/ 


2 2 2 
saa  = (sse  + ssn  + ssh  ) 


1/2 


sa  = 


sse  — ssn  — , ssh  r- 

e + n + h 

ssa  ssa  ssa 


= saa  e + sab  n + sac  h . 

The  discussion  for  camera  b is  precisely  analogous, 
yielding  the  desired  unit  vector 

sb  = sba  e + sbb  n + sbc  h . 

In  determining  the  solution  point  the  unit  vectors  sa 
and  sb  in  the  ray  directions  must  be  known  at  the  same  instant. 
Since  solutions  are  desired  at  equally  spaced  time  points,  a 
method  of  interpolating  between  two  camera  rays  is  needed. 

Although  there  are  timing  flash  marks  on  each  roll  of 
film  which  allow  the  time  at  which  each  frame  was  exposed  to 
be  determined  very  closely,  the  shutters  of  the  two  cameras  are 
not  synchronized.  Moreover,  the  picture-taking  rates  of  the 
two  cameras  are  not  quite  equal. 

Assume  that  we  have  two  successive  frames  from  either 
one  of  the  cameras  taken  at  times  tas  and  tbs  where  tas  < tbs. 
Let  the  film  frame  coordinates  of  the  moving  point’s  image  be 
(ppa,  qqa)  at  time  tas  and  (ppb,  qqb)  at  time  tbs. 

Figure  44  shows  vectors  a and  b to  the  image  points  of 
the  moving  object  point  at  times  tas  and  tbs  respectively.  Now 
to  get  what  would  have  been  the  ray  direction  to  the  moving 
point’s  image  at  time  ts  (tas  < ts  < tbs) , we  assume  that  the 
ray  turns  with  constant  angular  velocity  through  the  angle  th  in 
time  interval  tbs  - tas.  That  is,  the  angle  ths  turned  through 
in  time  ts  - tas  is  ths  = ■—  th.  See  Figure  45. 
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Figure  44. 


b 


Figure  45. 

Figure  45  shows  vectors  a,  b,  and  the  angles  th  and  ths. 

Vector  c is  that  vector  which  picks  out  on  the  film  frame  the 

image  point  position  that  would  have  occurred  had  a picture  been 

taken  at  time  ts.  (The  assumptions  that  the  moving  point  remains 

in  the  plane  of  its  rays  a and  b from  time  tas  to  time  tbs,  and 

that  it  moves  in  such  a way  that  its  angular  velocity  with  respect 

to  the  camera  focal  point  is  constant  are  in  general  not  exactly 

ture,  but  if  the  film  frame  rate  i£  high  enough  to  faithfully 

record  the  experiment,  they  should  be  good  assumptions.  The 

simple  error  statistics  on  the  program  output  seem  to  indicate 

the  above  is  a good  interpolation.) 

‘ 


1 
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Referring  to  Figure  44,  we  have 
a = f = ppa  r - qqa  J 
b = 1 + ppb  T + qqb  J 


f = /f/ 


, „2 

2 

2 

1/2 

a = /a/  = 

(f 

+ ppa 

+ qqa  ) 

2 

1/2 

b = /b/  = 

(f 

+ ppb 

+ qqb  ) 

th  = sin  * (/r  x s/) 

ths  = -ts  ~-tas — th 
tbs  - tas 

tht  = th  - ths . 

Let  rs  be  a unit  vector  in  the  direction  of  c 
Figure  46. 


b 


Figure  46 


So 


rs  = A r + B s where  A and  B must  be  found. 
Dotting  through  by  r gives : 

r.rs  = Ar.r  + Br.s 


or  cos  ths  = A + B cos  th. 

Dotting  through  by  s gives : 

ri".  s-Ar.s  + Bs.s 
or  cos  tht  = A cos  th  + B 


cos 

ths 

cos  thi 

cos 

tht 

1 1 

1 1 

cos  th  | 

| cos 

th 

1 | 

cos 

(ths)  - cos  th  cos  tht 

1 

2T7 
- cos  th 

cos 

(th 

- tht)  - cos  th  cos 

sin2  th 

sin 

th  sin  tht 

. 2 
sm 

th 

sin 

tht 

sin 

th 

1 1 

cos  ths i 

Icos 

th 

cos  tht 1 

. 2 
sm 

th 

cos  tht  - cos  ths  . cos  th 

T^~T 

sm  th 
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sin  ths 
sin  th 


rs  = A r + B s 


— sin  tht  — , sin  ths  — 

rs  = — > ZXL — r + I rrr  s 

sin  th  sin  th 


where  rs  is  a unit  vector  in  the  direction  of  the  image  of  the 
moving  point  at  time  ts. 

We  need  only  scale  rs  so  that  its  tip  just  touches  the 
film  frame.  See  Figure  47. 


(p»q) 


Figure  47. 


The  length  of  f is  f , i.e.,  f = /f/ 
c = xk  rs  where  xk  must  be  found. 


xk  rs  . 7 = 


f (the  component  of  c in  the  direction  of 
the  camera  optical  axis  must  equal  the 
focal  length,  f,  of  the  camera. 


xk  « 


rs  . ~E 
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and 


c = xk  rs. 

The  coordinates  of  the  image  point  at  time  ts  are: 

(p,q)  where  p = c . I 
and  q = c . J . 

This  essentially  concludes  the  discussion  of  the  math 
modeling  used  in  the  program.  The  rest  of  the  program  is  principally 
concerned  with  bookkeeping. 

It  is  worth  mentioning  that  the  ray  interpolation  probably 
could  have  been  done  more  economically.  However,  an  additional 
complication  with  the  runs  is  that  some  slight  camera  movement  is 
experienced  during  periods  of  peak  sled  acceleration.  It  is  very 
slight  on  the  great  majority  of  runs.  There  are,  however,  one  or 
two  runs  in  which  camera  vibration  was  rather  severe.  In  order 
to  deal  with  this  problem,  additional  analysis  (not  reported  in 
detail)  has  been  done  to  allow  frame-by-frame  reorienting  of  the 
cameras.  To  insert  this  correction  into  the  program  would  require 
a different  ray  interpolation  routine. 

The  assumption  of  an  absolutely  rigid  camera  mounting 
appears  to  be  justified  in  view  of  the  data  obtained  so  far.  But 
reorienting  of  the  cameras  frame-by-frame  might  give  some 
improvement.  Also,  taking  the  solutions  from  the  program  as  it 
presently  stands  and  actually  going  back  and  computing  the  angular 
accelerations  of  the  moving  points  at  each  camera  would  provide 
an  interesting  check  on  the  ray  interpolation  assumptions  presently 
in  the  program.  It  might  even  allow  a reinterpolation  and  a better 
final  solution. 

In  short,  there  are  a number  of  areas  that  remain  to  be 
investigated  in  this  problem,  perhaps  with  profit,  perhaps  not. 
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3.2.1  Program  SLED 


Program  SLED  uses  four  subroutines.  Its  block  diagram 
appears  in  Figure  48. 


Figure  48.  Block  Diagram  of  Program  SLED. 

A brief  overview  of  the  program  operation  will  be  given 
before  the  more  detailed  discussion.  The  main  program  SLED  reads 
in  the  input  data,  scales  it  where  necessary,  and  writes  it  out. 

In  broad  outline  this  input  data  consists  of  location  and 
orientation  data  for  each  of  the  two  cameras  and  the  position 
coordinates  of  certain  reference  points  rigidly  affixed  to  the 
sled.  Also,  the  film  frame  coordinates  of  the  moving  points  are 
input  for  each  of  the  cameras  for  each  film  frame.  The  film  frame 
coordinates  of  some  of  the  fixed  reference  points  are  also  input, 
but  at  present  no  use  is  being  made  of  this  information  in  program 
SLED.  The  time  increment  desired  between  solution  points  is 
input  and  a run  identification  number.  The  exact  reference  times 
for  at  least  two  frames  are  read  in  for  each  camera.  Subroutine 
TMSGN  (Time  Assign)  then, by  linear  interpolation,  assigns  to 
each  of  the  other  film  frames  its  correct  time. 


Subroutine  ALL  is  then  called.  It  computes  the  time, 
ts,  at  which  each  of  the  four  moving  points  is  to  be  determined. 
It  then  finds  that  pair  of  film  frames  from  each  camera  whose 
times  bracket  the  present  ts  value.  Then,  for  the  first  targeted 
point  ALL  calls  FIND  once  for  each  camera  to  get  the  interpolated 
rays  toward  that  point  at  time  ts.  ALL  then  calls  SOLVE  to  find 
the  position  of  the  point  at  time  ts.  This  is  done  in  succession 
for  each  of  the  four  moving  points  on  the  test  subjects.  Then, 
ts  is  increased  by  the  desired  time  increment  between  solution 
points  and  the  whole  process  repeated  again.  This  continues  until 
a time,  ts,  is  requested  which  exceeds  the  time  of  the  last  frame 
input  for  one  of  the  cameras.  Throughout  this  process  the  sum  of 
the  two  ray  solution  miss  distances  are  stored  for  each  of  the 
four  moving  points.  Also  stored  is  the  sum  of  the  squares  of 
these  miss  distances.  After  the  four  moving  points  have  been 
located  for  the  last  value  of  ts,  the  sample  mean  and  sample 
standard  deviation  of  the  miss  distances  is  computed  for  each  of 
the  four  points. 

The  principal  symbols  used  in  the  program  will  now  be 
defined.  This  first  list  consists  of  those  variables  appearing 
in  COMMON.  COMMON  consists  solely  of  camera  location  and 
orientation  data. 


RAA  s 

, 

RAB 

RAC  > 

[ Camera  a 

position  coordinates  (feet). 

FA 

Camera  a 

focal  length  (feet). 

THA 

Camera  a 

optical  axis  azimuth  (radians) . 

PHA 

Camera  a 

optical  axis  elevation  (radians) . 

GA 

Camera  a 

tilt  angle  in  radians.  Not 

presently 

used. 

THIA 

Azimuth  of  the  camera  a film  frame  i factor 
(radians) . 

PHIA 

Elevation 
(radians) . 

of  the  camera  a film  frame  i vector 
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i 


1 


Camera  b positon  coordinates  (feet) . 


I 


RBA 

RBB 

RBC 


I 


FB 

Camera 

b 

focal  length  (feet). 

THB 

Camera 

b 

optical  axis  azimuth  (radians) . 

PHB 

Camera 

b 

optical  axis  elevation  (radians) 

GB 

Camera 

used. 

b 

tilt  angle  in  radians.  Not  pres 

THIB 

Camera 

b 

i axis  azimuth  (radians) . 

PHIB 

Camera 

b 

r axis  elevation  (radians) . 

The  next  part  of  the  symbol  list  consists  principally  of 
the  dimensioned  variables  used  by  the  program.  A short  discussion 
follows  this  list. 


NA 

NFRA(I)  , 
1=1,  NA 

NB 

NFRB(I)  , 
1=1,  NB 

TA(I)  , 

I = 1,  NA 

TB(I)  , 

I = 1,  NB 

XA(I,J) , 

I = 1,  NA; 
J = 1,  8 

YA  (I,J) , 
1=1,  NA; 
J = 1,  8 

XB  (I,J), 
1=1,  NB; 
J = 1,  8 

YB  (I,J) , 

I = 1,  NB; 
J = 1,  8 


Number  of  frames  of  camera  a data. 

Frame  number  of  the  ith  frame  of  camera 
a data. 

Number  of  frames  of  camera  b data. 

Frame  number  of  the  ith  frame  of  camera 
b data. 

Time  at  which  camera  a shot  ith  frame 
(seconds) . 

Time  at  which  camera  b shot  ith  frame 
(seconds) . 

The  horizontal  film  frame  coordinate  of 
the  jth  image  point  on  the  ith  film  frame 
of  camera  a (counts:  1000/counts/inch). 

The  vertical  film  frame  coordinate  of 
the  jth  image  point  on  the  ith  film  frame 
of  camera  a (counts*  1000/counts/inch). 


The  analogous  quantities  for  camera  b. 
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Up  to  10  camera  a frame  numbers  and 
their  associated  times  may  be  read  in 
and  used  by  TMSGN  to  assign  the  correct 
times  to  all  the  camera  a film  frames 
(TREFA  (I)  values  in  seconds). 


The  analogous  quantities  for  camera  b. 
In  practice  only  two  reference  frame 
numbers  and  times  need  be  used. 

The  sample  mean  and  the  sample  standard 
deviation  of  the  two-ray-solution  miss 
distances  for  each  of  the  four  targeted 
points  (inches) . 


The  position  coordinates  of  up  to  10 
reference  points  rigidly  attached  to 
the  sled  (feet). 

The  above  list  can  be  further  clarified  by  discussing 
those  quantities  in  the  above  list  that  pertain  only  to  camera  a. 

The  cameras  are  brought  up  to  speed  shortly  before  the  sled  begins 
to  accelerate  in  order  that  the  picture-taking  rate  of  each  camera 
be  constant  during  the  actual  experiment.  Thus,  the  first  frame 
of  camera  a containing  data  of  interest  might  be  frame  80. 

Assume  there  are  33  frames  of  camera  a data  numbered  consecutively 
from  80  to  112 . 

Then  NA  = 33 

NFRA(l)  = 80 

NFRA ( 2 ) = 81 


NFRA(33)  = 112 

Now  there  are  eight  points  on  each  of  these  frames  whose 
coordinates  are  read  on  the  film  frame  reader.  The  program  stores 
these  as  XA  (I,J),  YA(I,J)  1=1,  33;  J = 1,  8.  Thus,  the  fifth 
point  read  on  the  ninth  frame  would  have  coordinates  XA  (9,5), 

YA  (9,5).  This  is  on  frame  number  88,  i.e.,  NFRA  (9)  = 88. 


NREFA(I) , 

I = 1,  10 

TREFA  (I) , 
I = 1,  10 

NREFB(I) , 

I = 1,  10 

TREFB (I) , 

I = 1,  10 

SMEAN  (J) , 
J = 5,  8 

SSDEV(J)  , 

J = 5,  8 

XX(10) 

YY (10) 

ZZ (10) 
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The  cameras  do  run  at  a constant  rate  once  they  have  come 
up  to  speed.  There  are  timing  flashes  of  light  superimposed  on 
the  negative  from  which  the  time  an  individual  frame  was  shot  can 
be  determined  quite  accurately.  Suppose  that  frame  80  is  found 
to  have  been  taken  at  time  0.142  seconds  and  frame  112  at  0.208 
seconds.  Then  the  program  would  store 

NREFA(l)  = 80  NREFA  (2)  = 112 

TREFA(l)  = 0.142  TREFA  (2)  = 0.208. 

Subroutine  TMSGN  then  assigns  the  correct  time  to  each 
frame  of  camera  a data  by  linear  interpolation  and  stores  them 
in  TA (I) ,1=1,  33,  i.e., 

TA  (1)  = 0.142 

TA  (2)  = 0.144 


TA  (32)  = 0.206 

TA  (33)  = 0.208. 

It  is  not  necessary  to  use  the  first  and  last  frame  of 
data.  Had  frames  88  and  99  been  used  as  reference  frames,  TMSGN 
would  linearly  extrapolate  to  those  frames  outside  the  interval. 

Precisely  the  same  discussion  applies  to  those  variables 
associated  with  camera  b.  Since  the  two  cameras  are  not  synchronized, 
the  only  connection  between  their  data  is  the  time  flash  which  is 
inscribed  on  both  film  strips  simultaneously.  Thus,  camera  a and 
camera  b in  general  have  different  numbers  of  frames  and  in 
general,  no  two  frames  from  a and  b are  snapped  simultaneously. 
Moreover,  it  is  generally  true  that  a test  time  ts  at  which 
the  positon  of  the  moving  points  are  to  be  output  is  not  coincident 
with  that  of  any  frame  from  either  camera. 

The  input  for  the  program  will  now  be  discussed. 


The  first  four  cards  contain  the  location  and  orientation 


information  for  camera  a and  b as  follows: 

Cards  1 and  2:  RAA,  RAB,  RAC,  FA,  THA,  PHA,  GA 

THIA,  PHIA 

Format  (4F10 . 0/5F10 . 0) 

Cards  3 and  4:  RBA,  RBB , RBC , FB,  THB , PHB , GB,  THIB, 

PHIB 

Format  (4F  10 . 0/5F10 . 0) . 

The  input  of  camera  a and  camera  b film  frame  data  is 
terminated  by  signal  cards  at  the  ends  of  their  respective  data 
decks.  The  numbers  NA  and  NB  are  determined  during  the  read- in 
of  this  data.  The  data  for  each  film  frame  requires  two  cards. 

The  first  card  for  each  frame  contains  the  frame  number  and  then 
the  film  frame  coordinates  for  the  images  of  the  four  fixed  points. 
The  second  card  contains  the  same  frame  number  and  then  the  film 
frame  coordinates  for  the  four  moving  points.  NA  is  incremented 
by  one  as  the  two  cards  for  each  new  frame  of  camera  a data  are 
read  in. 

Camera  a data  input  for  each  frame 

NFRA(NA)  , ( (XA(NA,I)  , YA(NA,I)),I  = 1,  8) 

Format  (15,  8F7.0/5X,  8F7.0). 

Termination  cards  for  camera  a input  consists  of  two 
cards,  the  first  with  9's  in  columns  1-5,  the  second  blank. 

99999 

A BLANK  CARD 

Camera  b data  input  for  each  frame  is 

NFRB(NB),  ( ( XB  (NB , I ) , YB  (NB  , I)  ) , I = 1,  8) 

Format  (15,  8F7.0  / 5X,  8F7.0). 

Termination  is  the  same  as  camera  a. 

99999 

A BLANK  CARD. 
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Next  the  posi  tlon  coordinates  of  some  reference  points  are 
read  in.  NHP , the  number  of  reference  points,  is  defined  as  they 
are  read  in.  Again,  the  read-in  is  terminated  by  a signal  card. 

The  format  is 

XX (NHP) , YY (NHP) , ZZ(NHP) 

Format  (3F10.0). 

The  termination  card  has  100000.0  punched  in  the  first 
ten  columns  and  is  otherwise  blank. 

The  reference  times  and  the  reference  frame  numbers  for 
camera  a are  then  input  to  the  program.  NTMA,  the  number  of 
reference  times  for  camera  a,  is  determined  when  the  read-in  is 
terminated  by  a card  with  99999  in  columns  12-16. 

Camera  a data  input  for  each  reference  time  and 
frame  number  is 

TREFA  (NTMA) , NREFA  (NTMA) 

Format  (F10.0,  16) 

The  termination  card  is 

99999  (in  columns  12-16)  . 

The  camera  b time  reference  input  is  precisely  analogous 
and  is  terminated  in  the  same  way. 

TREFB  (NTMB) , NREFB  (NTMB) 

Format  (F10.0,  16) 

with  termination  card 

99999  ( in  columns  12-16) . 

The  last  input  card  contains  the  run  identification  number 
DESIG  and  the  time  increment  DT  (in  seconds)  desired  between  position 
solutions  for  the  moving  points.  The  format  is 

DT,  DESIG 

Format  (2F10.0). 

Figure  49  indicates  the  various  sections  of  a typical 
input  deck . 
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DT  DESIG 


99999 


Time 


TREFB  (NTMB)WREFB (NTMB) 


99999 


TREFA  (NTMA)  NREFA  (NTMA) 

100000.0 


/ Reference  Point  Deck 


XX (NHP)  YY  (NHP)  ZZ  (NHP) 


99999 
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) , 


99999 

/!fe®ra  A Film  Frame  Data 
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THB  PHB  GB  THIB  PHIB 
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THA  PHA  GA  THIA  PHIA 


RAA  RAB  RAC  FA 


u 


Figure  49.  Typical  Data  Deck  Input  to  Program  SLED 


mr- 

There  is  a great  deal  of  output  presently  taken  from 
the  program.  Much  of  this  is  debug  output  which  can  now  be 
eliminated  and  this  will  be  ignored.  Rather,  that  part  of  the 
input  that  is  read  out  will  be  discussed  as  well  as  those  sections 
of  output  pertaining  directly  to  the  displacement  versus  time 
solutions  of  the  four  moving  points  on  the  dummies.  Also,  the 
sample  mean  and  sample  standard  deviation  checks  on  the  goodness 
of  the  solutions  will  be  discussed. 

When  this  program  was  in  its  inception,  it  was  assumed 
that  there  would  be  either  points  read  from  each  film  frame  for 
each  camera,  the  first  four  being  the  image  points  of  four  refer- 
ence points  rigidly  attached  to  the  sled  and  the  second  four  being 
the  four  targeted  points  on  the  dummies.  It  was  planned  to  run 
the  two-ray  solution  for  all  eight  points,  hoping  that  any  bias 
in  the  orientation  of  the  cameras  would  be  revealed  by  the  solutions 
for  the  four  fixed  points.  Unfortunately,  when  some  fixed  points 
are  visible  from  one  camera,  they  are  often  obscured  from  the 
other.  For  that  reason,  the  sets  of  four  fixed-points  that  can  be 
read  from  each  camera  throughout  the  experiment  differ  from  one 
another. 

However,  although  the  two  camera  solutions  cannot  be  run 
for  these  differing  sets  of  four  fixed  points,  this  data  would 
still  be  necessary  if  a frame-by-frame  reorientation  of  the  cameras 
should  be  attempted  in  the  data  reduction . In  the  two  runs  where 
camera  vibration  was  severe,  this  approach  should  be  tried.  In 
all  the  other  cases  in  which  camera  vibration  is  slight,  the 
solutions  obtained  have  been  quite  good  under  the  present  program 
assumption  that  there  is  no  camera  vibration.  What  improvement 
could  be  achieved  in  these  cases  by  the  incorporation  of  a frame- 
by-frame  reorientation  scheme  into  the  program  is  simply  not  known 
at  this  time. 

As  the  program  presently  stands,  it  makes  no  use  whatso- 
ever of  the  first  four  image  points  read  in  for  each  frame  of  each 
camera.  Nor  is  any  use  made  of  position  coordinates  of  the 
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reference  points  read  in.  This  information  is  being  retained  with 
a view  as  to  modifications  yet  to  be  made  to  program  SLED.  For 
the  above  reasons,  the  sample  mean  and  sample  standard  deviation 
read-out  for  the  first  four  points  are  recorded  as  0.0  since  no 
solutions  are  attempted  for  these  points. 

Figure  50  is  a tab  page  showing  some  of  the  output. 

The  first  two  lines  are  the  camera  location  and  orientation  data 
for  cameras  a and  b.  Position  coordinates  and  focal  lengths 
are  in  feet  and  all  the  angles  in  radians. 

The  bracketed  data  immediately  below  consists  of  the 
position  coordinates  of  each  of  the  reference  points  input.  These 
are  in  feet  and  are  given  with  respect  to  the  program  right-hand 
coordinate  system  as  are  the  camera  positions  above.  The  reference 
times  are  in  seconds.  Part  of  the  frame  data  for  camera  a is 
output.  There  are  three  lines  for  each  frame  of  data.  Taking  the 
first  frame  as  an  example,  the  "1"  identifies  this  as  the  first 
frame  of  camera  a data,  the  "80"  is  the  film  frame  number,  and 
the  .141864E  + 00  is  the  time  in  seconds  at  which  it  was  snapped. 
The  next  line  gives  the  film  frame  image  coordinates  of  the  four 
fixed  points  read  for  camera  a and  the  following  line  contains 
the  coordinates  of  the  four  targeted  points  on  the  dummies.  All 
these  film  frame  image  coordinates  are  output  in  feet.  After  all 
the  film  frame  data  for  both  cameras  a and  b has  been  output, 
the  time  increment,  dt,  desired  between  solution  positions  for 
each  of  the  four  points  is  printed  out  in  seconds.  Refer  to  the 
complete  listing  beginning  on  page  F-15  to  see  this. 

Figure  51  displays  that  page  of  the  tab  on  which  the  sample  mean 
and  the  sample  standard  deviation  for  each  of  the  four  targeted  points  (point 
numbers  5 through  8)  are  provided.  These  quantities  are  the  mean  and  the 
standard  deviation  of  all  the  miss  distances,  d,  in  the  two-ray  solutions  for 
each  of  the  moving  points.  They  are  output  in  inches. 
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The  last  section  of  output  gives  the  positions  of  the 
moving  points  on  the  subjects  versus  time.  See  Figure  52. 

The  position  versus  time  data  is  first  output  for  points 
5 and  6,  the  forehead  and  chin  of  the  subject  in  the  passenger 
seat.  The  first  quantity  is  the  time  in  seconds.  The  next  three 
numbers  are  the  x,  y,  and  z coordinates  of  point  five  in  inches. 
The  next  three  numbers  are  the  coordinates  of  point  six  in 
inches.  The  coordinates  of  the  moving  points  are  output  with 
respect  to  the  original  left-hand  experimental  coordinate  system 
affixed  to  the  sled.  The  run  number  is  just  an  identifier.  Then 
points  5 and  6 are  specified.  All  this  data  is  also  punched  out 
on  cards  for  input  to  other  programs.  The  last  three  numbers  of 
each  line  of  this  output  are  not  punched  on  the  cards. 

These  numbers  are  debug  output.  The  first  is  the  distance 
between  the  point  pair  5 and  6 to  the  left.  The  second  and  third 
are  the  miss  distance,  d,  that  occurred  in  the  computation  of  the 
positions  of  points  5 and  6. 

A complete  listing  of  program  SLED  is  presented  in 
Appendix  F. 
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SECTION  4 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  processes  developed  and  applied  during  the  period 
of  performance  have  satisfied  the  basic  requirements  for  which 
they  were  developed.  There  is  still  much  that  could  be  accomplished 
to  refine  the  programs  for  the  purposes  of  improving  cost 
effectiveness  and  user  confidence  in  the  solutions  generated. 
Continuing  review  is  recommended  for  the  purpose  of  accomplishing 
the  following. 

A.  Modify  program  HIFPD  to  include  a computerized 
analysis  of  differences  between  the  raw  and  smoothed  solution 
points . 

B.  Modify  program  HIFPD  to  provide  computerized  cal- 
culations of  conversion  factors  for  the  various  anthropometric 
points  being  tracked. 

C.  Investigate  the  desirability  of  employing  smoothing 
techniques  other  than  the  moving  arc  quadratic  fit  currently 
being  used. 

D.  Determine  accuracy  requirements  and  criteria  for 
analyzing  the  accuracy  of  solutions  prior  to  applications  of  these 
processes  to  future  studies. 
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DESCRIPTION  OF  PROGRAM  INPUT  DATA  AND  PARAMETER  CODES 


I. 


PROGRAM  SETUP  CARDS 

A)  The  first  card  in  the  setup  deck  must  contain  the  date  in 
columns  1 to  10;  for  example,  12  FEB  74  or  FEB  11,74  (only 
one  date  card  per  job) . 

B)  The  following  four  cards  are  required  for  each  test  in 
the  computer  job: 


Card  Number  1 
Column  Format 

1-5  A5 

6 11 

7 11 

8 11 

9 11 


11  11 


13  11 


14-16  13 

17-19  13 


Data  Description 
Test  number 

IRX — flag  controlling  polarity  of  X-axis 
data-  blank  or 

0  no  change 

1  change  sign  of  X-axis  data 

IPR — flag  controlling  input  data  and 
difference  printout  -blank  or 

0  print  data 

1  omit  printout 

I TYPE — flag  controlling  input  data-blank 

or  0 read  and  process  all  8 variables 

1 read  and  process  only  variables 

1,2,7,  and  8 

IPL — flag  controlling  linear  velocity  and 
and  acceleration  data  -blank  or 

0  print  and  plot  data 

1  print  only 

2  omit  all  linear  data 

IPA — flag  controlling  angular  velocity 
and  acceleration  data  (shoulder-hip  and 
head  pt.l,  head  pt.2) — blank  or 

0  print  and  plot  data 

1  print  only 

2  omit  all  angular  data 

IPC — flag  controlling  variable  displacement 
with  respect  to  sled  data  -blank  or 

0  print  and  plot  data 

1  print  only 

2  omit  these  computations 

First  frame  included  in  displacement 
with  respect  to  the  sled  plot  (if  blank, 
the  first  available  frame  is  plotted) . 

Last  frame  in  displacement  with  respect 
to  the  sled  plot  (if  blank,  the  last 
frame  is  plotted) . 
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Card  Number  2 (Continued) 


Column 

Format 

Data  Description 

20-21 

12 

The  number  of  sets  (M)  of  linear  velocity  and 
acceleration  to  be  computed. 

23-24 

211 

Variable  code  of  variable  and  reference  respectively 
(see  variable  codes  in  Item  II)  for  first  set  of 
linear  data;  for  example, 

21 sled  relative  to  range  or 

32 hip  relative  to  sled 

26-27 

211 

Same  as  above  for  the  2nd  set  of  linear  data 

29-30 

211 

Same  as  above  for  3rd  set 

Repeat 

motion 

the  above 
(Format : 

for  each  of  the  M (maximum  of  12)  sets  of  linear 

IX,  211  for  each  set) . 

Card  Number  3 

1-10 

F10.0 

Time  calibration number  of  frames  per  second. 

May  be  left  blank  if  film  speed  is  500  frames 
per  second . 

11-20 

F10.0 

SLED  calibration  in  counts  per  foot 

21-30 

F10.0 

HIP  calibration  in  counts  per  foot* 

31-40 

F10.0 

KNEE  calibration  in  counts  per  foot* 

41-50 

F10.0 

SHOULDER  calibration  in  counts  per  foot* 

51-60 

F10.0 

ELBOW  calibration  in  counts  per  foot* 

61-70 

F10.0 

HEAD  POINT  1 calibration  in  counts  per  foot 

71-80 

F10.0 

HEAD  POINT  2 calibration  in  counts  per  foot 

NOTE: 

The  decimal  must  be  punched  in  the  above  data  fields  unless 
the  data  are  integer  and  are  right  justified. 

Card  Number  4 

1 

1L 

9 in  column  1 to  indicate  the  end  of  test  input 

NOTE:  Cards  1,  2,  and  3 are  placed  in  front  of  the  test  deck 

and  card  4 is  placed  after  the  last  frame  in  the  test. 


C)  The  last  card  in  the  input  deck  (before  the  end  of  job  card) 
contains  the  word  "END"  in  columns  1 to  3. 


II. 


VARIABLE  CODE  IDENTIFICATION 


The  following  code  versus  variable  name  list  is  used  through- 
out the  program  and  card  2 in  Item  1-B: 


Code 

1 

2 

3 

4 

5 

6 

7 

8 


Name 

Range 

Sled 

Hip 

Knee 

Shoulder 

Elbow 

Head  Point  1 
Head  Point  2 


III. 


CARD  FORMATS  FOR  THE  TEST  INPUT  DATA  CARDS  FOR  ITYPE=0 
Card  Number  1 


Column 

Format 

Data  Description 

1 

11 

Card  identification  (must  be  1) . 

2-5 

14 

Frame  number 

6-12 

17 

X reading  in  counts  for  Range  data 

13-19 

17 

Z reading  in  counts  for  Range  data 

20-26 

17 

X for  Sled 

27-33 

17 

Z for  Sled 

34-40 

17 

X for  Hip 

41-4  7 

17 

Z for  Hip 

48-54 

17 

X for  Knee 

55-61 

17 

Z for  Knee 

Card  Number  2 

1 11 

Card  identification  (must  be  2) 

2-5 

14 

Frame  Number 

6-12 

17 

X reading  in  counts  for  Shoulder  data 

13-19 

17 

Z reading  in  couts  for  Shoulder  data 

20-26 

17 

X for  Elbow 

27-33 

17 

Z for  Elbow 

34-40 

17 

X for  Head  Point  1 

41-47 

17 

Z for  Head  Point  1 

48-54 

17 

X for  Head  Point  2 

55-61 

17 

Z for  Head  Point  2 
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IV.  CARD  FORMATS  FOR  THE  TEST  INPUT  DATA  CARDS  FOR  ITYPE=1 


Card  Number  1 


Column 

Format 

Data  Description 

1 

11 

Card  : 

identi f ication 

(must  be  1)  . 

2-5 

14 

Frame 

number 

6-12 

17 

X 

reading  in  counts 

for  Range  data 

13-19 

17 

Z 

reading  in  counts 

for  Range  data 

20-26 

17 

X 

for 

Sled 

27-33 

17 

Z 

for 

Sled 

34-40 

17 

X 

for 

Head  Point  1 

41-47 

17 

Z 

for 

Head  Point  1 

48-54 

17 

X 

for 

Head  Point  2 

55-61 

17 

z 

for 

Head  Point  2 

NOTE: 

For  I TYPE 

= 1 

. , only  1 data  card 

is  read  for  each 

frame . 


V.  GENERAL  COMMENTS 

A)  If  there  are  any  errors  in  frame  or  card  identification 
numbers,  error  statements  will  be  printed  at  the  top 

of  the  first  output  page  for  the  test  and  all  com- 
putations after  the  listing  of  the  input  data  will  be 
deleted. 

B)  A maximum  of  150  frames  (MAXN)  will  be  read  for  each 
test.  If  the  test  input  deck  contains  more  than  150 
frames,  only  the  first  150  will  be  processed.  This 
could  be  changed  by  changing  MAXN  and  the  array 
dimensions  in  the  program. 

C)  If  the  calibration  factor  for  a variable  is  missing, 
flag  I CAL ( J)  is  set  equal  to  zero  and  that  variable 
will  be  deleted  from  the  analysis. 

D)  An  eleven  point  quadratic  least-square  fit  is  used 
throughout  the  program.  This  could  be  changed  by 
changing  the  value  of  "NP"  in  the  program. 

E)  The  CALCOMP  plot  abscissa  and  ordinate  scales  for  the 
velocity  and  acceleration  data  are  determined  by  the 
plot  subroutine  (CPLT) . The  X (abscissa)  and  Z 
(ordinate)  scales  for  the  variable  relative  to  the  sled 
data  are  constants.  The  X-scale  ranges  from  -1.4  to 
2.6  feet  for  IRX  = 0 and  -3.6  to  0.4  feet  for  IRX  = 1 
and  Z scale  ranges  from  0.0  to  4.0  feet.  Any  displace- 
ments outside  this  range  will  be  set  equal  to  the 
limiting  value. 
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HYGE  Program  Setup  Cards 


1)  DATA  card,  used  once  per  job:  | | | | | — | — ) — ) — | — | 

Col . 1 col .10 

e.g.  12  FEB  74  (first  card  in  input  deck). 

2)  The  following  four  cards  are  required  for  each  test: 

I CARD  No . 1 / TITLE  Card  --80  columns  of  alphanumeric 

information. 


Col  .1 


80 


I CARD  No I 2~| 


TEST  No. 


FIRST  FRAME 
2 


IRX  Q IPR  Q I TYPE  □ I PL  □ IPA  □ IPC  ["1 

6 7 8 9 11  13 


LAST  FRAME 
5 6 7 


M Sets  -*■ 


oo  ca  m m m m m rn 

2627  29  30  3233  3536  3739  4142  4445  4748 


50 


10 

m 

0 51 


11 


53  54  56  5 


12 

L '/  ^ (Maximum  of  12  sets;  usually  3 or  4.) 


CARD  NO.  3 | Calibration  Data .( frames/second  or  counts/foot) 


TIME: 


1 10 
SHOULDER*  | 


SLED  | | HIP  * f I KNEE*  | | 

11  20  21  30  31  40 


ELBOW* 


41  50 

* Must  be  0 or  blank  for  "I TYPE  = 1" . 


i I HEAD  PT.l  I I HEAD  I I 

51  60  61  70  PT.2  ft  80 


3) 


CARD  NO . 4~~) 

Card  No.  4 has  a "9"  in  column  1 (insert  after  the  last  frame 

in  the  data  deck). 

END  Card,  used  once  per  job:  Punch  "END"  in  columns  1 to  3. 
(last  card  in  input  deck;  before  "END  OF  JOB"  Card) . 
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V. 


General  Comments: 


A)  If  there  are  any  errors  in  frame  or  card  identification  numbers, 
error  statements  will  be  printed  at  the  top  of  the  first  output  page 
for  the  test  and  all  computations  after  the  listing  of  the  input  data 
will  be  deleted. 

B)  A maximum  of  150  frames  (MAXN)  will  be  read  for  each  test.  If  the 
test  input  deck  contains  more  than  150  frames,  only  the  first  150 
will  be  processed.  This  could  be  changed  by  changing  MAXN  and 
the  array  dimensions  in  the  program. 

C)  If  the  calibration  factor  for  a variable  is  missing,  flag  ICAL(J)  is 
set  equal  to  zero  and  that  variable  will  be  deleted  from  the  analysis. 

D)  An  eleven  point  quadratic  least  square  fit  is  used  throughout  the 
program.  This  could  be  changed  by  changing  the  value  of  "NP"  in 
the  program. 

E)  The  CALCOMP  plot  abscissa  and  ordinate  scales  for  the  velocity 
and  acceleration  data  are  determined  by  the  plot  subroutine  (CPLT). 
The  X (abscissa)  and  Z (ordinate)  scales  for  the  variable  relative 

to  the  sled  data  are  constants.  The  X scale  ranges  from  -1.4  to 
2.  6 feet  for  IRX  = 0 and  -3.  6 to  0.  4 feet  for  IRX  = 1 and  the  Z scale 
ranges  from  0.  0 to  4.  0 feet.  Any  displacements  outside  this  range 
will  be  set  equal  to  the  limiting  value. 
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appendix  b 

PREPARATION  OF  DATA  FOR  INPUT  TO  HIFPD 
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1 

STANDARD  PRACTICE  PROCEDURE 
93291-02-1  (Rev.  "I",  4 Oct.  ’76) 

This  procedure  describes  the  process  to  be  employed  in  the 
reduction  of  photometric  data  acquired  during  experimental  tests 
conducted  on  the  Horizontal  Impulse  Accelerator  during  the 
Restraint  Systems  Dynamics  (RSD)  investigation. 

1 . Source  Documents 

1.1  16  mm  motion  picture  film  from  cameras  mounted: 

a.  Onboard  Lateral  view 

b.  Onboard  Oblique  view 

c.  Offborad  Lateral  view  (Backup  for  "A") 

d.  Offborad  Downtrack  view 

e.  Offboard  Overhead  view 

1.2  Restraint  System  Dynamic  Pretest  Anthropometric 
Measurements  data  sheet. 

1.3  Sled  Coordinate  System  and  Camera  Scheme,  (RSD). 

2 . Equipment  Required 

2.1  Producers  Service  Corp.  film  digitizer  coupled  to 
Teletype  w/tape  punch. 

2.2  Teletype  Terminal  w/tape  reader  and  telephone 
coupler. 

2.3  IBM  Keypunch,  Type  029. 

3.  Editing  Requirements 

The  following  sequence  of  operations  shall  be  followed  to 
determine  if  all  necessary  data  have  been  adequately 
recorded. 

3.1  Preview  Onboard  lateral  view  film  to  determine  if 
the  following  points  are  observable  throughout  the 
test  period  from  first  motion  of  the  subject  with 
respect  to  the  seat  through  return  to  rest  of  the 
subject  with  respect  to  the  seat.  (Refer  to 
Figure  B-l) . 
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a. 


Forward  fiducial,  LH  side  of  seat 


b.  Aft  fiducial,  LH  side  of  seat 

c.  Fiducial  over  hip  (or  fiducial  over  iliac  creat) 

d.  Fiducial  on  knee 

e.  Fiducial  on  shoulder 

f.  Fiducial  on  elbow 

g.  Fiducial  on  jaw  hinge  (or  head  mount) 

h.  Accelerometer  at  origin  of  9 TAP  accelerometer 
assembly. 


3.2 


Determine  time  - film  frame  table.  Assume  time  to 
be  0.000  seconds  at  upper  edge  of  frame  in  which 
the  synchronizing  flash  is  first  observed,  set 
frame  counter  to  00000.  Interpolate  time  lapse 
between  this  point  and  next  0.1-second  time  pulse. 
(Refer  to  Figure  B-2)  . Set  this  value  to  te^. 

Using  the  same  method,  determine  the  time  lapse 
between  the  last  .01  second  time  pulse  prior  to  the 
frame  at  which  the  subject  returns  to  rest  with 
respect  to  the  seat  (Frame  "n"  and  the  start 
(or  top)  of  this  frame.  Let  this  value  be  t^. 

Count  the  0.01  second  time  intervals  between  the 
first  pulse  after  frame  000  and  the  last  pulse 
prior  to  frame  "n".  Multiply  the  number  of 
intervals  by  0.01  seconds  and  let  this  value  be 
tc  (calibrated  time) . 

The  total  elapsed  time  can  now  be  calculated  by 
summing  t ^ + t _ + tc. 


117 


Figure  B-2.  Time  Lapse  Determination. 


The  film  speed  (or  frames  per  second) , which  is 
required  as  input  to  the  computer  program  "HIFPD" 
can  now  be  calculated: 

Film  Speed  = Frame  Count  at  Frame 

el  e e2 


3.3  Determine  Conversion  Constants  using  the  "widths" 
information  recorded  on  the  "Restraint  Systems 
Dynamics  Pretest  Anthropometric  Measurements"  data 
sheet,  (Figure  B-3)  the  "Sled  Coordinate  System 
and  Camera  Scheme"  (Figure  B-4)  and  the  PSC 
digitizer. 


Read  the  coordinate  of  the  fiducials  mounted  on  the 
side  of  the  seat.  Record  these  on  the  Calibration 
Form  (Figure  B-5)  . 


Calculate  resultant  scalar  dimension  between  the 
fiducials  on  the  panel  (s  .)  and  the  seat  (s  .): 

c pi  si 


Spi  = (Xpl  ' Xp2>2  + (*pl  - *p2)2 


Ssi  = lxsl  - xs2>  2 + '^1  ‘ ^s2 1 2 ' 

Calculate  conversion  constants  for  these  dimensions : 


s . 

f = Hi 

p 11.96875 

12 


In  like  manner: 


f 


s 


12.0 

12 


Where  f and  f 
P s 

of  counts/foot. 


are  conversion  constants  in  terms 
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RESTRAINT  SYSTEM  DYNAMICS 
Photometric  Calibration  Form 


TEST 


Fiducial  Readings: 

Seat  Fwd : x1 

Seat  Aft.  x2 

Panel  Fwd. 

Panel  Aft.  x„ 


Resultant  Scalar  Dimensions: 


s 


SI 


counts 


yl 

y2 

yl 

y2 


Pi 


counts 


Conversion  Factors: 
s_ 


si 


s 


ct/ft 


fP  = 


Focal  Point  to  Seat  Distance: 

24.062 


11.96875 

12 


Ct/ft. 


s 


-(>) 


24.062 


Conversion  Constant  Input  File: 
film  speed 


hip 


knee 


tshoulder 


elbow 


jaw 


accel. 


ss 

sr 

+ 8 - 

s 

ss 

i 

+ 8 - 

s 

Ss 

I 

+ 8 - 

s 

s 

s 

+ 8- 

s 

ss 

s 

+ 8- 

s 

Ss 

x f 


x f 


X f = 

s 


X f 


X f = 

s 


X f 


frames/second 

counts/foot 

counts/foot 

counts/foot 

counts/foot 

counts/foot 

counts/foot 

counts/foot 

TEST 


Figure  B-5.  Photometric  Calibration  Form. 
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The  distance  from  the  focal  point  of  the  camera  to 
the  vertical  plane  of  the  edge  of  the  seat  can  now 
be  calculated  (Reference  Figures  B-6  and  B-7) . 

Since  the  distance  from  the  transparent  panel  to  the 
seat  edge  has  been  measured  to  be  24.062-inches 
and  the  approximate  distance  from  the  focal  plane 
of  the  camera  to  the  seat  edge  plane  was  measured 
to  be  60.625-inches  then  the  distance  from  the  focal 
point  of  the  lens  to  the  lexan  panel,  the  seat 
edge,  and  thus  the  seat  can  be  calculated  by 
similar  trianges  (Reference  Figure  B-7)  . 


s . x 
si 


s . x 
si 


Knowing  the  measured  distance  sg  - sp  to  be 

24.062  inches,  thus  s = s - 24.062,  by 

P s 

substituting  we  have: 


s - 24.062 
s 


and 


1 - 24.062 


f 

fs 


f 

f 


JL 

By  subtracting  from  both  terms  and  adding 


24.062 


s 


s 

to  both  terms : 


, _ 24.062 

fs  = ss 
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Seat 


T 


Fiducials 
Mounted  to 
this  Side 


24-1/16 


Figure  B-6.  Schematic  of  Camera-Reference  Fiducials. 


1 


and  by  multiplying  both  terms  by  sg  and  by 

1 


we  have  that: 


Since  we  assume  the  focal  length  of  the  camera 
lens  to  be  precisely  as  stated  by  the  manufacturer, 
the  final  calculated  s could  not  be  expected  to 

b 

coincide  exactly  with  the  measured  distance. 
However,  the  calculated  distance,  sg,  would  be 
expected  to  be  in  close  agreement  with  the 
measured  distance. 

Once  the  camera  to  seat  fiducial  distance  has  been 
determined,  conversion  factors  can  be  calculated 
for  the  planes  of  motion  of  the  antropometric 
points  to  be  tracked. 

If  we  assume  that  any  point,  p,  will  travel  in  a 
plane,  y = pt  then: 


conversion  factor  in  plane  y = pt 
( ft/count) 

conversion  factor  in  plane  of  seat 
f iducials 

distance  along  optical  axis  from 
focal  point  to  the  plane  y = pt. 


where : f 


Pt 


’Pt 


24.062 

-&■) 
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and: 


sg  = distance  along  optical  axis  from  focal 
point  to  the  plane  of  the  seat  fiducals 

S . is  determined  by: 
pt 

s = s„  + 8 inches 
pt  s 

where  w„  is  the  width  measurement  between  the 
P 

anthropometric  points  on  opposite  sides  of  the 
subject  as  noted  on  the  pretest  anthropometric 
measurements  form,  and  s =8  inches  = Sa  and 
the  subject  is  assumed  to  be  seated  such  that  his 
plane  of  symmetry  coincides  with  the  seat  (£,  and 
that  the  antropometric  points  on  the  left  and 
right  sides  of  the  subject  are  symmetrical  about 
that  plane. 

All  of  the  preceding  calculations  are  provided  for 
on  the  photometric  calibration  form. 

4 . Film  Digitizing 

The  anthropometric  points  are  tracked  by  digitizing  the 
"x"  and  "y"  coordinated  of  each  point  as  projected  from 
each  frame  during  the  impact  response  period.  The  film 
must  be  mounted  in  such  a manner  that  the  subject  appears 
to  face  the  left  as  viewed  on  the  projection  screen. 
Because  of  the  smoothing  routines  utilized  in  the 
computer  processing  of  the  data,  care  must  be  taken  to 
start  digitizing  the  data  at  least  fifteen  (15)  frames 
before  the  subject  starts  to  move  relative  to  the  seat  and 
it  must  be  continued  at  least  fifteen  (15)  frames  after 
the  subject  returns  to  rest  relative  to  the  seat. 


_E 

2 
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Data  points  are  digitized  as  follows: 


Line  1 : Frame  Number 

a.  Forward  seat  reference  fiducial;  x,y 

b.  Aft  seat  reference  fiducial;  x,y 

c.  Hip  (or  iliac  crest)  fiducial;  x,y 

d.  Knee  fiducial;  x,y 

Line  2 : Frame  Number 

e.  Should  fiducial;  x,y 

f.  Elbow  fiducial;  x,y 

g.  Jaw  (or  accel.  pack  mount)  fiducial;  x,y 

h.  9 TAP  accelerometer  reference;  x,y 

The  above  data  are  to  be  listed  on  the  teletype  printer 
and  punched  into  paper  tape  on  the  teletype.  Care  must 
be  taken  to  punch  a series  of  rubouts  before  and  after 
digitizing  the  data. 


APPENDIX  C 
PROGRAM  HIFPD 
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APPENDIX  F 

PROGRAM  SLED  AND  TYPICAL 
LISTING  OF  SOLUTIONS 
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ERRATA  STATEMENT 

AMRL-TR-78-94 , "Photometric  Methods  for  the  Analysis  of  Human 
Kinematic  Responses  to  Impact  Environments",  October  1978  is 
revised  as  follows: 

1.  Page  11  - First  complete  paragraph,  sixth  and  seventh 
sentences, 

"The  resultant  displacement  . . . and  in  G's.", 
is  changed  to  read, 

A moving  eleven  point  quadratic  least  square  fit  is  then 
applied  to  these  smoothed  x and  z-axis  displacement  data  to  obtain 
the  x and  z components  of  velocity.  Next  this  same  smoothing 
routine  is  applied  to  these  x and  z-axis  velocity  data  to  compute 
the  x and  z components  of  acceleration.  The  resultant  displace- 
ment, velocity,  and  acceleration  data  are  then  computed  using 
these  smoothed  x and  z component  data. 

2.  Page  15  - The  equation  used  to  calculate  resultant  dis- 
placement is  changed  to  read: 

RES  CD  = Ad  CD  2 + ZD  (I)  2 

The  second  sentence  following  this  equation, 
"Subroutine  DERIV1  ...  in  each  case.",  is  changed  to 
read. 

Next  subroutine  DERIV1  is  called  to  compute  the  X and  Z com- 
ponents of  velocity  from  the  X and  Z displacement  data  and  then 
to  compute  the  X and  Z components  of  acceleration  from  the  X and 
Z velocity  data.  In  each  case  an  eleven  point  quadratic  smoothing 
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function  is  applied  to  the  input  displacement  or  velocity  data.  ! 

The  velocity  components  are  stored  in  arrays  VXCI)  and  VZ(I)  and 
the  acceleration  components  in  arrays  AX  CD  and  AZ CD . The 
resultant  velocity  and  acceleration  data  are  computed  as  follows: 

VEL  CD  = / VXCD2  + VZ  CD  2 

ACC  Cl ) * /ax  CD  2 + AZ(I}2 

3.  Appendix  C - as  follows: 

Page  130  

Insert  line  125: 

2,  VX  (302)  , VZC302),  AX(302),  AZC302) 

Insert  line  185: 

2,  CXXCD1)  , VXCD)  , (XX  (1,2),  AX(1)),  CZZ(1,1), 
VZCl)),  CZZ(.1,2),  AZ  Cl)  ) 

Page  143  

Change  Line  4480  to: 

CALL  DERIVl  CTf  XD,  VX,  N,  NP,  1) 

Change  Line  4500  to: 

CALL  DERIVl  CT,  VX,  AX,  N,  NP,  2) 

Insert  Line  4485: 

CALL  DERIVl  (T,  ZD,  VZ,  N,  NP , 1) 

Insert  Line  4505: 

CALL  DERIVl  (T,  VZ , AZ , N,  NP,  2) 

Insert  Line  4525: 

VEL ( I ) = SQRT (VX (I)  * VX(I)  + VZ(I)  *VZ(I)) 

Insert  Line  4527: 

ACC  (I ) = SQRT  (AX  (I)  * AX(I)  + AZ(I)  * AZ  (I)  ) 

4.  Appendix  D - Disregard  tabular  listings  and  plots  of 

Resultant  Velocity  and  Resultant 
Accelerations . 


